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Abstract 

This paper extends the Gotay-Nester and the Dirac theories of constrained 
systems in order to deal with Dirac dynamical systems in the integrable case. In- 
tegrablc Dirac dynamical systems arc viewed as constrained systems where the 
constraint submanifolds are foliated, the case considered in Gotay-Nester theory 
being the particular case where the foliation has only one leaf. A Constraint Al- 
gorithm for Dirac dynamical systems (CAD), which extends the Gotay-Ncstcr 
algorithm, is developed. Evolution equations are written using a Dirac bracket 
adapted to the foliations and an abridged total energy which coincides with the 
total Hamiltonian in the particular case considered by Dirac. The interesting 
example of LC circuits is developed in detail. The paper emphasizes the point 
of view that Dirac and Gotay-Nester theories are dual and that using a com- 
bination of results from both theories may have advantages in dealing with a 
given example, rather than using systematically one or the other. 

^Corresponding author. Fax: -1-54-221-424-5875 
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1 Introduction 

The main purpose of this paper is to study an extension of the Dirac theory of 
constraints (Dirac [1950, 1958, 1964]) and the Gotay-Nester theory (Gotay, Nester, 
and Hinds [1978]; Gotay and Nester [1979]) for the case in which the equation of 
motion is a Dirac dynamical system on a manifold M, 



where D C TM ®T*M is a given integrable Dirac structure on M, which gives a 
foliation of M by presymplectic leaves, and £^ is a given function on M, called the 
energy. One of the features of systems like (1.1) is that they can represent some 
constrained systems where the constraints appear even if the Lagrangian is nonsin- 
gular, as it happens in the case of integrable nonholonomic systems. We will work 
under explicit hypotheses about regularity, locality and integrability, although some 
of our results can be applied in more general indicated in the paragraph 

Hypotheses, below in this section. 

In the Gotay-Nester theory the starting point is an equation of the form 



where w is a presymplectic form on M. In fact, in the Gotay-Nester theory the more 
general case where d£{x) is replaced by an arbitrary closed 1-form a is considered, 
but this will not be relevant here. 

Example: Euler-Lagrange and Hamilton's equations. Let L: TQ — )• M be 

a Lagrangian, degenerate or not. Let £{q, v,p) = pv — L{q, v) and let u E Q^{TQ © 
T*Q) be the presymplectic form co = dq^ A dpi on the Pontryagin bundle M = 
TQ®T*Q. Then the Euler-Lagrange equations are written equivalently in the form 
of equation (1.2) with x = {q,v,p). In fact, we have, in local coordinates. 



(x,x) © d£{x) G D. 



(1.1) 



uj{x){x, ) = d£{x) 



(1.2) 



Hq,v,p)d(f A dpi 



q'dpi - pidq"- 



(1.3) 



d£ 




(1.4) 



(1.5) 



Using equations (1.3)^(1.5) we can easily see that (1.2) becomes 



q' = v' 




which is clearly equivalent to the Euler-Lagrange equations. The case £{q,v,p) = 
H{q,p), where H is a given Hamiltonian on T*Q of course gives Hamilton's equa- 
tions. 



1 Introduction 3 



The idea of using the Pontryagin bundle to write such fundamental equations of 
physics appears in Livens [1919]; Skinner [1983]; Skinner and Rusk [1983a, b]; Cen- 
dra, Marsden, Pekarsky, and Ratiu [2003b]; Yoshimura and Marsden [2006a, b]. 

Equation (1.2) is equivalent to equation (1.1) in the case in which the Dirac 
structure on M is the one naturally associated to co, denoted Di^, in which case 
the foliation of M has obviously only one presymplectic leaf. In general, equation 
(1.1) may be considered as a collection of equations of the type (1.2), one for each 
presymplectic leaf of the Dirac structure. However, in order to study the initial 
condition problem this approach might not be appropriate, because different initial 
conditions might belong to different presymplectic leaves and therefore correspond 
to different equations of the type (1.2), which is not the usual situation. 

The algorithm of the Gotay-Nester theory generates a finite sequence of sec- 
ondary constraint submanifolds M D Mi ^ . . . D Mc- The final constraint subman- 
ifold Mc has the property that every solution curve of (1.2) is contained in Mc (in 
fact, it is the smallest submanifold having that property). Equations of motion are 
given by restricting the variable {x,x) in equation (1.2) to TMc, and existence of 
solution curves x{t) for a given initial condition xq G Mc is guaranteed, under the 
condition that the kernel of 

uj\x)\T,Mc: T,Mc^T*M 

has locally constant dimension. An important point in the Gotay-Nester approach 
is that the equations defining the secondary constraint submanifolds are written in 
terms of uj (see formula (2.5)), which makes the whole algorithm invariant under a 
group of transformations preserving uj. 

In order to solve (1.1), we will develop a constraint algorithm that we call CAD 
{constraint algorithm for Dirac dynamical systems). The CAD extends the Gotay- 
Nester algorithm and gives a sequence of secondary foliated constraint submanifolds. 
The secondary foliated constraint submanifolds can be written in terms of the Dirac 
structure using formula (5.6), which generalizes (2.5). 

Let M C be an embedding of M into a symplectic manifold {N,fl), in such 
a way that the presymplectic leaves of M are presymplectic submanifolds of A^ and 
let, by slight abuse of notation, £ : N ^ M denote an arbitrary extension of the 
given energy on M. This kind of framework can be naturally constructed in many 
examples. Then one can consider, at least locally, that M is a primary foliated con- 
straint submanifold defined by equating certain constraints to 0, which generalizes 
the situation of the Dirac theory of constraints. One of the results of this paper 
establishes that there exists a Dirac bracket whose symplectic leaves are adapted 
to the foliations of the primary and final foliated constraint submanifolds. We also 
prove that the equations of motion on the final foliated constraint submanifold can 
be nicely written in Hamiltonian form with respect to this Dirac bracket and the 
abridged total energy, which we introduce in section 3.4. 

Our extension of the Gotay-Nester and Dirac theories has especially simple and 
interesting features when applied to LC circuits. We show that the algorithm stops 
either at Mi or M3, and concrete physically meaningful formulas for the first three 
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constraint submanifolds and for the evolution equations can be given for this case. 
Other geometric treatments of LC circuits can be found, for instance, in Bfoch and 
Crouch [1999]; Yoshimura and Marsden [2006a, b]; Cendra, Marsden, Ratiu, and 
Yoshimura [2011]. 

Systems of the type (1.1) are important also in the not necessarily integrable 
case since they represent, for instance, the Lagrange-d'Alembert equations of non- 
holonomic mechanics. More references will be given in section 4.1. We remark that 
even though the main focus of this paper is the case where the Dirac structure 
is integrable, many results, most notably the constraint algorithm CAD and the 
treatment of nonholonomic systems in section 5, are proven for a not necessarily 
integrable Dirac structure. 

Gotay-Nester's theory is more geometric than Dirac's and provides a certain 
notion of duality or correspondence between the two, in the spirit of the dualities 
between submanifolds and the equations defining them, manifolds and their ring 
of functions, Hamiltonian vector fields and their Hamiltonian or Poisson brackets 
and the collection of their symplectic leaves. Emphasizing this point of view in a 
common geometric language and showing that a combination of results from both 
theories may have advantages in a given example, rather than using one or the other, 
is an aspect of the paper. 

From this point of view, the notion of second class constraints is related to the 
notion of second class constraint submanifold, the latter being determined in this pa- 
per as a submanifold that is tangent to a second class subbundle of a certain tangent 
bundle. They are symplectic submanifolds (Sniatycki [1974]; Bursztyn [2010]). The 
first class constraint submanifolds are coisotropic submanifolds (Bursztyn [2010]). 

The presence of primary and secondary first class constraints, which in Dirac's 
work is connected to the principle of determinacy and to the notion of physical 
variables, is explicitly related here to certain manifolds M and M in a commuta- 
tive diagram, in Theorem 3.19. The diagonal arrow in this diagram, which is a 
submersion, is dual to the diagonal arrow in the diagram in Lemma 3.17, which is 
a surjective Poisson algebra homomorphism. In some examples, instead of apply- 
ing the Dirac method to write equations of motion, it may be simpler to proceed 
in two stages. First apply the Gotay-Nester algorithm to find the final constraint 
submanifold, which, as proven in the Gotay-Nester theory, coincides with the sub- 
manifold defined by equating the final constraints obtained in the Dirac algorithm 
to 0. Then switch to the Dirac viewpoint and write equations of motion using the 
Dirac bracket. Also, in order to calculate the Dirac bracket at points of the final 
constraint submanifold, using the symplectic leaf that contains that submanifold 
may be easier than using Dirac's expression. 

The point of view of IDEs. In order to study the problem of finding solution 
curves of general Dirac dynamical systems (1.1) the more general theory of Implicit 
Dijferential Equations may be useful. An Implicit Differential Equation (IDE) on a 
manifold M is written as 




(1.6) 
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A solution of (1.6) at a point x is a vector (x, v) € T^M satisfying (1.6). A solution 
curve, say x{t), t £ {a,b), must satisfy, by definition, that {x{t),x{t)) is a solution 
at x{t) for all t G (a, b). 

Basic questions such as existence, uniqueness or extension of solutions for a given 
initial condition are not completely answered yet, although many partial results 
have been established for certain classes of IDE. One approach is to use a constraint 
algorithm, which consists of a decreasing sequence of constraint submanifolds M ^ 
Ml 3 . • ■ 3 Mc defined as follows, 

Mk+i := {x G Mk I there exists {x,v) £ TMk such that ip{x,v) = 0}, 

with Mq := M. This algorithm, which obviously uses only the differentiable struc- 
ture regardless of any other structure which may be present, like presymplectic, 
Poisson or Dirac structures, represents a differential geometric aspect underlying 
the algorithms of Gotay-Nester, Dirac or CAD. To ensure that each is a sub- 
manifold and that the algorithm stops after a finite number of steps one may choose 
to assume certain conditions of the type "locally constant rank" conditions. Then 
the original IDE is reduced to an equivalent ODE depending on parameters on the 
final constraint submanifold Mc- In fact, by construction, is characterized by 
the property that it is the smallest submanifold that contains all solutions curves of 
the given IDE. Therefore, if is empty, there are no solution curves. 

Dirac's original work has a wide perspective from the physical point of view, 
with connections to classical, quantum and relativistic mechanics. However, from 
the point of view of abstract IDE theory, and in a very concise way, we may say that 
a combination of the Dirac and the Gotay-Nester methods shows how to transform a 
given Gotay-Nester equation (1.2) into an equivalent ODE depending on parameters 
(3.24) on a final constraint submanifold, while what we show in this paper is how 
to transform a given Dirac dynamical system (1.1) into an equivalent ODE (7.11) 
depending on parameters on a final foliated constraint submanifold. 

Some more comments on the connection between Dirac's and Gotay-Nester ideas 
and IDE are in order. We can compare Dirac [1950]; Gotay, Nester, and Hinds [1978]; 
Gotay and Nester [1979]; Rabier and Rheinboldt [1994]; Pritchard [2003]; Cendra 
and Etchechoury [2006] , to see how the idea of applying a constraint algorithm works 
in different contexts. In Cendra and Etchechoury [2006], one works in the realm of 
subanalytic sets; in Gotay, Nester, and Hinds [1978] and Gotay and Nester [1979] 
one works with presymplectic manifolds; in Pritchard [2003] one works with complex 
algebraic manifolds; Dirac [1950] uses Poisson brackets; in Rabier and Rheinboldt 
[1994] some degree of differentiability of the basic data is assumed, and, besides, 
some constant rank hypotheses are added, essentially to ensure applicability of a 
certain constant rank theorem. Some relevant references for general IDEs connected 
to physics or control theory, which show a diversity of geometric or analytic methods 
or a combination of both are Carihena and Rahada [1984]; Gracia and Pons [1991, 
1992]; Mendella, Marmo, and Tulczyjew [1995]; de Leon and de Diego [1996]; Ibort, 
de Leon, Marmo, and Martin de Diego [1996]; Marmo, Mendella, and Tulczyjew 
[1997]; Delgado-Tellez and Ibort [2003]. 
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Hypotheses. As for regularity, we will work in the C°° category. Throughout 
the paper we assume that the subsets appearing in sequences of the type 
M D Ml 5 • • • 2 Mc that are generated by some constraint algorithm, are sub- 
manifolds regularly defined by equating some functions (constraints) to 0, and that 
the sequence stops. More regularity conditions like Assumptions 3.8, 3.18, 3.26, 
3.27, 7.1, and some others will be introduced when needed along the paper. 

Our results will be of a local character, but for some of them, like the notion of 
second class submanifolds, it is indicated how to define them globally. The usage of 
local coordinates is almost entirely avoided and basic facts in symplectic geometry 
or Poisson algebra arguments are used instead. 

The condition of integrability of the Dirac structure appearing in (1.1) has its 
own interest. However, the CAD does not assumes integrability and can be applied 
for instance to general nonholonomic systems. On the other hand, in the non- 
integrable case certain brackets that we study would not satisfy the Jacobi identity, 
and will not be studied in this paper. 

Structure of the paper. The first part of the paper, which includes sections 2 
and 3, contains a review of the Dirac and Gotay-Nester methods. 

In sections 4 to 7 we develop our main results, extending the Gotay-Nester and 
Dirac theories. Section 4 is devoted to a review of basic facts on Dirac structures 
and Dirac dynamical systems (Cendra, Marsden, Ratiu, and Yoshimura [2011]). The 
notion of Dirac structure (Courant and Weinstein [1988]; Courant [1990]; Bursztyn 
and Crainic [2009]) gives a new possibility of understanding and extending the the- 
ory of constraints, which is the main purpose of the present paper. In section 5 
we develop our CAD algorithm. We do it for the general case of a not necessarily 
integrable Dirac structure so that one can apply CAD to general noholonomic sys- 
tems. In section 6 we study the example of integrable nonholonomic systems, and 
LC circuits are viewed as a particular case. In section 7 we show how to extend the 
Dirac theory for the case of a Dirac dynamical system (1.1). 

2 A brief review of the Dirac and Gotay-Nester theories 

In this section we review, without proof, basic facts belonging to the Dirac and the 
Gotay-Nester theories. 

The Dirac method starts with a given submanifold M, called the primary con- 
straint submanifold, of a symplectic manifold (iV, Q) defined by equating the pri- 
mary constraints to 0, and a given energy £ : N ^ M.. In the original work of Dirac 
N = T*Q, is the canonical symplectic form, M is the image of the Legendre 
transformation and £ is the Hamiltonian defined by Dirac. However, locally, any 
symplectic manifold is an open subset of a cotangent bundle, therefore the origi- 
nal Dirac formalism can be applied to the, seemingly more general, case described 
above, at least for the local theory. 

An interesting variant of the original situation considered by Dirac is the fol- 
lowing. Consider the canonical symplectic manifold = T*TQ with the canonical 
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symplectic form 0,, and let the primary constraint be M = TQ © T*Q, canoni- 
cally embedded in via the map given in local coordinates {q,v,p,i') of by 
if{q,v,p) = {q,v,p,0). In particular, M is defined regularly by the equation u = 0. 
If Lv is the presymplectic form on M obtained by pulling back the canonical sym- 
plectic form of T*Q using the natural projection, then one can show that (p*0, = oj. 
The embedding ip is globally defined (see Appendix A for details). The number pv 
is a well-defined function on M and it can be naturally extended to a function on 
a chart with coordinates {q,v,p,u), but this does not define a global function on 
A^ consistently. On the other hand, it can be extended to a smooth function on A^ 
using partitions of unity and any such extension will give consistent equations of 
motion. In this paper we will not consider global aspects. In any case, for a given 
Lagrangian L : TQ — )• M we can take £ = pv — L{q, v). 

Remark 2.1. For a given presymplectic manifold (M, cj) one can always find an 
embedding ip into a symplectic manifold {P,Q) such that f*Q = oo. Moreover, this 
embedding can also be chosen such that it is coisotropic, meaning that '^{M) is 
a coisotropic submanifold of P (see Gotay and Sniatycki [1981/82]). However, we 
should mention that the embedding given above is not coisotropic. 

The Dirac and the Gotay-Nester algorithms can be studied independently. On 
the other hand, they are related as follows. For a given system (1.2) choose a 
symplectic manifold (A^, il) in such a way that (M, w) is a presymplectic submanifold 
and £ (using a slight abuse of notation) is an arbitrary extension of £^ : M — )• M to A^. 
Moreover, assume that M is defined regularly by a finite set of equations i^-'^'' = 0, 
i = 1 , . . . , ao , where each cjif^ is a primary constraint. The Dirac algorithm gives 

a sequence of secondary constraints 4>'l'\ i = l,...,afc, which defines regularly 
a sequence of secondary constraint submanifolds M^, k = l,...,c, by equations 

(k) 

= 0, i = l,...,a/j, which coincide with the ones given in the Gotay-Nester 
algorithm. 

2.1 A brief review of Dirac's theory 

Dirac's theory of constraints has been extensively studied from many different points 
of view and extended in several directions. Part of those developments in the spirit 
of geometric mechanics is contained in the following references, but the list is far 
from being complete, Sudarshan and Mukunda [1974]; Gotay, Nester, and Hinds 
[1978]; Gotay and Nester [1979, 1980]; Mukunda [1980, 1987]; Skinner [1983]; Skinner 
and Rusk [1983a, b]; Carihena and Rahada [1984, 1993, 1995]; Carihena, Gomis, 
Ibort, and Roman [1985]; Cantrijn, Carihena, Crampin, and Ibort [1986]; van der 
Schaft [1987]; Carihena, Lopez, and Roman- Roy [1988]; Carihena [1990]; Henneaux 
and Teitelboim [1992]; Krupkova [1994]; de Leon, de Diego, and Pitanga [1995]; 
Krupkova [1996]; Marmo, Mendella, and Tulczyjew [1997]; Ibort, de Leon, Marrero, 
and Martm de Diego [1999]; Barbero-Lihan and Muiioz-Lecanda [2009]; Crampin 
and Mestdag [2011a,b]. 

As we have explained above we will work in a general context of a given sym- 
plectic manifold {N,0,), M C A^ and £ : N — )■ M, where the primary constraint 
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submanifold M is regularly defined by equations (p- = 0, i = 1, . . . ,aQ on N. 

The Dirac constraint algorithm goes as follows. One defines the total energy £t = 
£+X'^Q-^(f)f'\ The preservation of the primary constraints is written {(f)f'\£T}{x) = 0, 
i = 1, . . . ,aQ, X ^ Mo, or 

{cl)f\£}{x) + Ajo){</.f \<Af }(x) = 0, i,j = 1, . . . ,ao, X e Mo. 

Then Mi is defined by the condition that x G Mi if and only if there exists A(o) — 
(A^Q^, . . . , A"qP such that the system of equations 4>f\x) = 0, {4>f'\£T}{x) = 0, 

i = l,...,ao, is satisfied. The submanifold Mi is defined by equations 0^-^^ = 0, 
i = 1, . . . , fli, where each is a secondary constraint, by definition. By proceeding 
iteratively one obtains a sequence Mo 5 Mi 5 • • • 5 M^- Then there are final 
constraints, say , i = 1, . . . , Cc, defining a submanifold Mc by equations cpl = 0, 
i = 1, . . . , ac, called the final constraint submanifold, and the following condition is 
satisfied: for each x G Mc there exists (A^g^, . . . , A^q^ such that 

{<pf\£}{x) + X\^.^{ckf\(t>f}{x) = 0, i = 1, . . . , ae, J = 1, . . . , ao. (2.1) 

For each x G Mc the space of solutions of the linear system of equations (2.1) 
in the unknowns A^^q^ is an affine subspace of M'"' , called PSx^ , whose dimension is a 

locally constant function d^'^\x) = oq — Tank{{(j)^-'^\ 0^.^^}(x)). One can locally choose 
d^^^ (x) unknowns as being free parameters and the rest will depend affinely on them. 
Then the solutions of (2.1) form an affine bundle A^^^^ over M^. After replacing 
A(o) G A^"^^ in the expression of the total energy, the corresponding Hamiltonian 
vector field, 

XeAx) = Xe{x) + A;o)X^(o)(x), (2.2) 

X G Mc, which will depend on the free unknowns, will be tangent to Mc- Its integral 
curves, for an arbitrary choice of a time dependence of the free unknowns, will 
be solutions of the equations of motion, which is the main property of the final 
constraint submanifold Ale from the point of view of classical mechanics. The lack 
of uniqueness of solution for a given initial condition in Mc, given by the presence 
of free parameters, indicates, according to Dirac, the nonphysical character of some 
of the variables. In our context the physical variables will be given a geometric 
meaning. 

Remark 2.2. Dirac introduces the notion of weak equality for functions on T*Q. 
Two such functions are weakly equal, denoted / ~ (7, if /|Mc = g\Mc- Then, for 
instance (p^^'^ w 0. If / ~ then / = v^cfJ'f' , for some functions on T*Q and 
conversely. Since we use the notion of a constraint submanifold, in particular the 
final constraint submanifold, we prefer not to use the notation ~. 

Now let us make some comments on the notions of first class and second class 
constraints. The rank of the skew-symmetric matrix ({(/)-'^\ (/)^'^^}(x)), j = l,...,ac, 
is necessarily even, say, 2s, and it is assumed to be constant on Mc, as part of 
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our regularity conditions; for our results we will assume later a stronger condition 
(Assumption 3.27). One can choose, among the 4>f \ i = l,...,ac, 2s functions 
Xj, j = 1, . . . , 2s , and replace the rest of the by t^l'^^ + a^Xj with appropriate 
functions in such a way that tpi = 0, Xj = 0) define Mc regularly and, besides, 
{il^i,il^i'}{x) = 0, {^pi,Xj}{x) = 0, det{{xj,Xj'}{x)) / 0, for = l,...,Oc- 
2s, = l,...,2s and x G Mc- The (^^'^^ are linear combinations with smooth 
coefficients of the Xj V'ij conversely. The functions Xj^ j = 1, . • . ,2s, are 
called second class constraints and the functions i = 1, . . . ,ac — 2s, are called 
first class constraints. 

More generally, any function p on T*Q satisfying p\Mc = 0, {p,il^i}\Mc = 0, 
Xj}|-^c = 0, is a first class constraint with respect to the submanifold Mc, by 
definition. Any function g on T*Q satisfying {(7, ■0j}|Mc = 0, {g.,Xj}\Mc = 0, is a 
first class function, by definition. For instance, the total energy £t is a first class 
function. 

Now define the energy £c in terms of V'i, Xji ^ = 1) • • • > Oc — 2s, j = 1, . . . , 2s, as 

£c = £ + X'iJi + p^Xj- 

The preservation of the constraints for the evolution generated by £c can be rewritten 
as {ipi,£c}{x) = 0, which is equivalent to {ilji,£}{x) = for all x G Mc, and 
{Xj,£c}ix) = 0, for all x G Mc- The latter is equivalent to 

{Xi,£}{x) + pHxi,Xj}{x) = 0, i,j = 1,. . . ,2s, 

for all X G Mc, which determines the p^ as well-defined functions on Mc- Then the 
solutions {p{x). A) form an affine bundle with base Mc and whose fiber, parametrized 
by the free parameters A, has dimension Oc — 2s. 

Any section {p{x),\{x)) of this bundle determines £c as a first class function. 
This means that X£^{x) G T^Mc, for each x G Mc, and therefore a solution curve 
of X^^ is contained in Mc provided that the initial condition belongs to Mc- The 
function £c is essentially the extended Hamiltonian defined by Dirac. 

Dirac defines an interesting bracket, now called the Dirac bracket, 

{F,Gr = {F,G} - {F,XiVnXj,G}, 

which is defined on an open set in T*Q containing Mc, where c*-^, which by definition 
is the inverse matrix of {Xi^Xj}: is defined. The Dirac bracket is a Poisson bracket 
and has the important property that for any function F on T*Q, the condition 
{F, Xj}* = 0, j = 1, . . . , 2s, is satisfied on a neighborhood of Mc, which implies 
that F = {F,£c} = {F,£c}*, for any function F- Besides, {V'j,V'i}* = 0, i,j = 
1, . . . , Cc — 2s, on Mc- Because of this, one may say that, with respect to the Dirac 
bracket, all the constraints Xj; j = 1) • • • > 2s and xpi, i = 1, . - - ,ac — 2s, behave like 
first class constraints, with respect to Mc- 

2.2 A brief review of the Gotay-Nester algorithm 

We will need just a few basic facts from the Gotay-Nester theory. 
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In order to find solution curves to (1.2) we can apply the general algorithm for 
IDEs described in the introduction, and the final IDE can be written 

ijj{x){x, ) = d£{x) 
{x,x) G TMc. 

Let ujc be the pullback of oj to Mc- If ooc is symplectic, one obtains the simpler 
equivalent equation 

ujc{x){x,) = d{8\Mc) (2.3) 

on Mc, which is in Hamiltonian form. However, one must be aware that if ojc is 
degenerate then (2.3) is not equivalent in general to (2.2). 
The equation on Mc, 

uj{x){X,) = d£{x), (2.4) 

where X G TMc, defines an affine distribution on Mc, more precisely, one has an 
affine bundle S^'^'^ with base Mc whose fiber si^^ at a given point x G Mc is, by 
definition, 

5^^) = {X e T^Mc I (2.4) is satisfied}. 

The equivalence between the Dirac and the Gotay-Nester algorithms can be 
made explicit as an isomorphism of affine bundles, as follows. The affine bundles 
S^""^ Mc and A^'^) Mc are isomorphic over the identity on Mc, more precisely, 
the isomorphism A*^"^) — )• S*-"^^ is given by A(o) — )■ X^,^, using equation (2.2). In 
particular, the rank of S^'^'^ is d^'^\ 

Describing the secondary constraints using w. The constraint manifolds M^ 
defined by the algorithm can be described by equations written in terms of the 
presymplectic form oj, which is a simple but important idea, because those equations 
will obviously be invariant under changes of coordinates preserving u. Depending 
on the nature of uj one may obtain analytic, smooth, linear, etc., equations, which 
may simplify matters in given examples. The context of reflexive Banach manifolds 
is used in Gotay, Nester, and Hinds [1978] and Gotay and Nester [1979]. 
The condition defining the subsets M^+i, /c = 0, 1, . . . namely, 

i„w(x) = d£{x) has at least one solution v G T^M^, 

is equivalent to 

d£{x) G {T^Mkt 

or, since {T.^M^f = {{T^MkYY , 

{d£{x),{T,MkT) = {{)}. (2.5) 

3 First and second class constraints, constraint sub- 
manifolds and equations of motion 

In this section we will describe some results on first class and second class constraints 
and constraint submanifolds and also equations of motion which will be useful for 
our extension of those notions, to be developed in section 7. 
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3.1 Linear symplectic geometry 

The following results about linear symplectic geometry are an essential part of many 
of the arguments that we use in the paper, since under our strong regularity assump- 
tions many of them are essentially of a linear character. 

Lemma 3.1. Let {E,n) be a symplectic vector space of dimension 2n, V E a 
given subspace. For a given basis ai, i = 1, . . . , r of V° , let X-i = a\, i = 1, . . . r. 
Then the rank of the matrix [ai{Xj)] is even, say 2s, and Xi, i = 1, . . . , r form a 
basis ofV^. Moreover, the basis ai, i = l,...,r can be chosen such that for all 

ai{Xj) = 6ij-.s, 1 < i < s 
ai{Xj) = -6i-s,j, s + 1 < i < 2s 
ai{Xj) = 0, 2s <i<r. 

Proof. Consider the subspace = {V°)^. By a well-known result there is a basis 
Xi, i = 1, . . . , r of such that for all j = 1, . . . , r 

0,{Xi,Xj) = 5ij-s, 1 <i < s 

Xj) = -6i-s,j, s + I < i <2s 
n{Xi,Xj) =0,2s<i<r 

then take = X^. The first part of the lemma is easy to prove using this. ■ 

Lemma 3.2. Let ai, i = 1, . . . , r be a basis of V° having the properties stated in 
Lemma 3.1. Then Xi, i = 2s + 1, . . . ,r form a basis ofVr\ . 

Proof. Let X = X'Xi be an arbitrary vector in . Now XXi G V CiV^ iS 
aj{X) = A*(aj(Xj)) = 0, j = l,...,r. Since the first 2s columns of the matrix 
[ai(Xj)] are linearly independent and the rest are zero, we must have A* = 0, for 
1 < i < 2s, and A*, i = 2s + 1, . . . ,r are arbitrary. This means that V H is 
generated by Xi, i = 2s + 1, . . . ,r. ■ 

Corollary 3.3. dimV" n = r - 2s. 

Proof. Immediate from Lemma 3.2. ■ 

Let oj be the pullback of 17 to y via the inclusion. Then {V, oj) is a presymplectic 
space. In what follows, the ^ and " operators are taken with respect to unless 
specified otherwise. 

Lemma 3.4. = Vr\V^. 

Proof. X G iff u{X,Y) = 0, Vy € F iff ^{X,Y) = 0, VF G This is 
equivalent to X gV r\V^. ■ 

Lemma 3.5. Let ji, i = 1, . . .r be a given basis of V° and let Yi = -yf, i = 1, . . . r. 

Let (3 €z E* be given. Then the following conditions are equivalent. 
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(i) /3(y'^) = 0. 
(a) The linear system 

/3{Yi) + X^j,{Yi) = (3.1) 
has solution A = (A"^ , . . . , A^) . 

Proof. Let us show that (3.1) has solution (A-*^, . . . , A^) iff the system 

(3{Xk) + fi'ai{Xk) = (3.2) 

has solution (^^, . . . , fi"^), where k,l = 1, . . . r and «/ is a basis satisfying the condi- 
tions of Lemma 3.1. Since Yi, i = 1, . . .r and X^, k = 1, . . . r are both bases of 
there is an invertible matrix [^4^ such that Xj^ = Al.Yi. Let [B^] be the inverse of 
[Al.], so Yi = B\Xi. Assume that (3.1) has solution A-', j = 1, . . .r. We can write 
(3.1) as 

P{Y,) + \m{Yj,Yi) = {),i = l,...,r. 
Using this we have that for k = 1, . . . , r 

= /3(4y,) + \m{Y,, A\Y,) = P{Xk) + \^niYj,Xk) 

= (3{Xk) + x'n{B'^Xi, Xk) = (3{Xk) + fi'n{Xi, Xk) 

where /x' = A-'iJj. This means that the system (3.2) has solution. The converse is 
analogous. Using this, Lemmas 3.2 and 3.4, and the form of the coefficient matrix 
[ai(Xfc)] in Lemma 3.1, the proof that (3.2) has solution /x = {fi^, . . . , fi^) iff (3{V^) = 
is easy and is left to the reader. ■ 

Lemma 3.6. Consider the hypotheses in Lemma 3.5. Then the solutions to 

ixu: = P\V (3.3) 

(if any) are precisely X = + X-'Yj, where (A^, . . . , A'') is a solution to (3.1). A 
solution to (3.3) exists if and only if j3{V'^) = 0. If u) is symplectic then (3.1) and 
(3.3) have a unique solution and if, in addition, (3^ G V , then A^ = 0, . . . , A'" = 
and (3^ coincides with X = (/3|y)^" defined by (3.3). 

Proof. Since Yj, j = 1, . . . ,r form a basis of we have that (A^, . . . , A'') is a 
solution to (3.1) iff {(3 + Xhj){y^) = iff /3 + X^-fj E iff /3« + X^Yj G V. Now, let 
X = (3'^ + X^Yj, where (A^, . . . , A*") satisfies (3.1). Then we have X G F as we have 
just seen and we also have 

ixuj = {ix^V = X'\V = (/3 + X^jj)\V = (3\V, 

since 7j, j = 1, . . . ,r generate V°. We have proven that X is a solution to (3.3). 
To prove that every solution X to (3.3) can be written as before, we can reverse 
the previous argument. Using this, it is clear that if oj is symplectic then (3.1) has 
unique solution, in particular, we have that det (7^(5^)) ^ 0. If, in addition, /S" G y 
then X^Yj = X — (3"^ € V. Since Y^-, j = 1, . . . , r is a basis of V^, using Lemma 3.4 
and the fact that = {0} we get that A^' = for j = 1, . . . , r. ■ 
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Corollary 3.7. Let A = {A | A satisfies (3.1)}. Then dim A = r — 2s = dimkero;. 

Proof, ker w = , which has dimension r — 2s from Corollary 3.3 and Lemma 3.4. 
On the other hand the dimension of the subspace of A satisfying (3.1) is clearly also 
r — 2s, since the coefficient matrix has rank 2s. ■ 

3.2 Poisson-algebraic and geometric study of constraints and con- 
straint submanifolds. 

We shall start with the constrained Hamiltonian system (A^, M), where (A^, Q) 
is a symplectic manifold, 8: N — )• M is the energy and Af C A^ is the primary 
constraint submanifold. The equation to be solved, according to the Gotay-Nester 
algorithm, is 

w(x)(X, ) = d£{x)\T^M, (3.4) 

where X G T^Mc and x £ Mc, Mc being the final constraint. Let ujc be the pullback 
of Vt via the inclusion of Mc in A^. Since Uc is presymplectic, kerwc is an involutive 
distribution. From now on we will assume the following. 

Assumption 3.8. The distribution kerwc has constant rank and defines a regular 
foliation Kc, that is, the natural map pxc- — ^ Mc, where Mc = Mc/Kc is a 
submersion. 

Lemma 3.9. The following assertions hold: 

(a) There is a uniquely defined symplectic form uJc on Mc such that p*j^^LJc = ujc- 

(b) Let X be a given vector field on Mc- Then there is a vector field X on Mc 
that is pKc-''^dated to X . 

(c) Let f £ J^{Mc) . Then there exists a vector field X on Mc such that X is pK^- 
related to Xj, and for any such vector field X the equality u}c{x){X, ) = d{p*j^^f){x) 
holds for all x G Mc . 

(d) Let XxQ G TxgMc. Then one can choose the function f € J-{Mc) and the 
vector field X in (c) in such a way that X(xq) = Xx^. 

Proof, (a) By definition, the leaves of the foliation Kc are connected submanifolds 
of Mc, that is, each p^^(z), z G Mc, is connected. For z G Mc, let x G Mc such that 
PKci^) = z. For A, B G TzMc, as pK^ is a submersion, there are A, B £ Mc such 
that TxPkA = TxPK^B = B. We define ujc{z){A,B) = uJc{x){A,B). To prove 
that this is a good definition observe first that it is a consistent definition for fixed 
X, which is easy to prove, using the fact that kera;c(x) = kerT-^pXc- Now choose a 
Darboux chart centered at x, say U xV , such that, in this chart, pK^ '■ U x V ^ U 
and iOc{x^,x'^) = Wc(x^), where uJc{x^,x'^) and Cddx^) are independent of {xi,X2). 
This shows that Oc is well defined on the chart. Using this and the fact that one 
can cover the connected submanifold p'^^^z) with charts as explained above, one can 
deduce by a simple argument that u)c{z) is well defined. 

(b) Let g he a Riemannian metric on Mc. Then for each x G Mc there is a 
uniquely determined X{x) G T^Mc such that X{x) is orthogonal to ker T-^pXc and 
TxPKc^{x) = ^"(2;), for all x G Mc. This defines a vector field X on Mc which is 
Pi^^ -related to X. 
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(c) Given / and using the result of (b) we see that there is a vector field X on 
Mc that is pi^^ -related to Xj. Then, for every x G Mc and every Yx G T^Mc, 

uJcix){X{x),Yx) = u:c{pkSx)) [Xj{pkSx)),TxPkJ^x) = df{pKSx)){TxPKcYx) 
= d{p*Kj){x){Yx). 

(d) One can proceed as in (b) and (c), choosing / such that {dj {pkSxq)))^ = 
TxqPKc^xo s-iid, besides, the metric g such that X^q is perpendicular to kerTxgPKc- 

■ 

Definition 3.10. (a) For any subspace A C J-{N) define the distribution Aa ^ TN 
by AA{x) = {Xf{x) \fGA}. 

(b) The space of first class functions is defined as 

= {/ G T{N) I Xf{x) G TxMc, for all x G Mc}. 

In other words, R^'^^ is the largest subset of J-{N) satisfying 

X G Mc. 

Remark 3.11. Dirac was interested in classical mechanics, where states are repre- 
sented by points in phase space, as well as in quantum mechanics where this is 
not the case. From the point of view of classical mechanics, among the constraint 
submanifolds and constraints the only ones that seem to play an important role are 
M, Mc and the constraints (i)f'\ (t>f^ defining them by equations 4>f^ = 0, cj^f' = 0, 
respectively. 

Lemma 3.12. (a) R^^'^ is a Poisson subalgebra of {J-{N),{ , }). 

(b ) Mc is an integral submanifold of A^(c) ■ Moreover, for any vector field X on 
Mc that is pKc-i"^l(^ted to a vector field Xj on Mc there exists a function f G R^'^^ 
such that f\Mc = p*p^ f and X = Xf\Mc. In particular, any vector field X on Mc 
satisfying X{x) G kera;c(x) for all x G Mc is px^-related to the vector field on 
the symplectic manifold Mc, which is associated to the function f = 0, therefore 
there exists a function f G R^'^\ which satisfies f\Mc = 0, such that X(x) = Xf(x), 
X G Mc. 

Proof, (a) Let f,g e R^^l Then Xf{x) and Xg{x) are both tangent to Mc at 
points X of Mc which implies that —X^f gy(x) = [Xf,Xg]{x) is also tangent to Mc 
at points of x of Mc. This shows that {/, g} G R^'^K It is easy to see that any linear 
combination of /, g and also fg belong to R^^'^ . 

(b) By definition A^(c) C TMc. We need to show the converse inclusion. Let 
XxQ G TxgMc, we need to find / G R^^^ such that Xf(xo) = Xxg. Choose the function 
/ and the vector field X on Mc as in Lemma 3.9, (d). Choose any extension of 
p*j^^f to a function g on A^. For each x G Mc, we can apply Lemmas 3.5 and 3.6 

with E := TxN, V := TxMc, (3 := dg{x), y,{x) := d(t)f\x), i = 1, . . . , r^. We obtain 
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that in a neighborhood [/ C of each point xq of Mc we can choose functions 
A*^^^(x), i = 1, . . . , Tc, such that 

X{x) = [dg{x) + Xl^{x)d(P^{x))\ 

for all X £ McClU. Let = c/(x) + Aj^)(x)(/>f ^(x), for ah j; E [/. Then we have 

that Xf^{x) = X{x), for aU x e McDU. 

Now, consider a partition of unity pi, i £ I, on N, where each pi is defined on 
an open set Ui, i G I. Let J C / be defined by the condition i G J if and only if 
Ui n Mc 7^ 0. Using standard techniques of partitions of unity and the above result 
one can assume without loss of generality that for each i G J there is a function fjj. 
defined on Ui such that Xj^ (x) = X{x), for all x G Mc H Ui. Let / = Ylii^j PifUi-, 
which can be naturally extended by on A^. Then it is easy to see, using the fact 
that fu^ix) = p*^ f{x) = g{x), for each x G Mc, that X{x) = Xf{x), for each 
X £ Mc, and in particular Xq = X{xo) = Xf{xo). ■ 

Lemma 3.13. (a) Each function f G R^'^^ is locally constant on the leaves of Kc 
therefore, since they are connected, for each f G i?^'^^ there is a uniquely determined 
f G F{Mc), called {pKc)*fj such that f\Mc = P*K^f- Moreover, the vector fields 
Xf(x), X G Mc, on Mc and Xj on Mc are pKc-'''e.lated. 

(h) For each f G T{Mc) there exists f G R'^'^^ such that f\Mc = P*k f "-nd the 
vector fields Xf{x), x G Mc, on Mc and Xj on Mc are px^-related. 

Proof, (a) Let / G we only need to show that / is constant on the leaves of 
Kc, which is equivalent to showing that df{x) \ keiujcix) = 0, for all x G Mc- For a 
given X G Mc, let Vx G keiudx); then using Lemma 3.12 (b) one sees that there is 
a function g G R^^'^ such that Vx = Xg{x). Then we have 

= iOc{x) {Xf{x), Xg{x)) = nix) {Xf{x),Xg{x)) = df{x)Xg{x). 

Now we shall prove that Xf and Xj are pi^-^-related. For each x G Mc and each 
Yx G TxMc, we have 

Uc{Xf{x),Yx) = n{Xf{x),Yx) = df{x){Yx) = d{p\J){x){Yx). 

Using this we obtain 

UJc{x){Xf{x),Yx) = Ujc{x){pK,{x)) {TxPK,Xf{x),TxPKYx) = df{pKA^)){TxPKYx), 

which shows that X j{pkX'^)) = TxPKcXf{x), because (Dc is symplectic and TxPRc^^x 
represents an arbitrary element of Tp^ (x)M. 

(b) To find / we choose a vector field X that is px^-related to according to 
Lemma 3.9 and then use Lemma 3.12 (b). ■ 

Definition 3.14. 

/W = {/ G I f\Mc = 0}, 

= {/ G R^"'^ \{f,h}e for all h G R^"'^}. 
Elements of I^^^ are called first class constraints. 
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Lemma 3.15. (a) /'■^^ is a Poisson ideal of R^^\ that is, it is an ideal of the ring 
R^''^ such that if f e 1^"^ then {f,h} e l(^), for all h E R^'^h 
(b) Zj(c)R^'^^ is a Poisson subalgebra of R^'^\ 

Proof, (a) Let f,g& I^'^^ and h E R^'^^ . Then it is immediate that f + g and kg 
belong to For any h E we have {/, h]\Mc = Xh{f)\Mc = 0. 

(b) Follows from (a), using basic Poisson algebra arguments. ■ 

Lemma 3.16. The following conditions are equivalent for a function f E i?'-'^^. 
fij /E%.)ii(^). 
(ii) f\Mc is locally constant. 
(Hi) Xf{x) E kera;c(x) for x E Mc- 

Proof. Assume (i). Then {/, h}\Mc = for all h E that is, df{x)Xh{x)\Mc = 
0. By Lemma 3.12, (b), we know that Xh{x) represents any vector in TxMc- We 
can conclude that f\Mc is locally constant, so (ii) holds. Now we will prove that (ii) 
implies (iii). Let f\Mc be locally constant. Then for all g E R^'^^ and all x E Mc, 

= Xg{f){x) = n{x){Xf, Xg)ix) = u;cix){Xf{x),Xg{x)). 

Since, again by Lemma 3.12, Xg{x) represents any element of TxMc, we can conclude 
that Xf{x) E kera;c(2^), so (iii) holds true. Now we will prove that (iii) implies (i). 
Assume that Xf{x) E keiujcix), x E Mc- Then for all g E i?^"^^ and all x E Mc, 

{g,f}{x) = n{x){Xg,Xf){x) = cJc{x){Xg,Xf){x) = 0, 

that is, {g, /} E Using this and the definitions, we see that f e Zj(,)R'^''\ U 

Lemma 3.17. The map (pKc)* '■ ^^^^ ~^ J'{Mc) defined in Lemma 3.13 is a surjec- 
tive Poisson map and its kernel is therefore there is a natural isomorphism of 
Poisson algebras {pkJ^Kc) : i^^^V^^'^^ ^ J'iMc). 

Proof. Surjectivity of (pKc)* the fact that its kernel is I^^^ follows immediately 
from Lemma 3.13 and the definitions. This implies that {pKc)*i{<^) is an algebra 
isomorphism. Also, using the definitions, for f,g € R^^^ and any x E Mc we can 
prove easily that 

{f,g}{x) = n{x){Xf{x),Xg{x))=ujc{x){Xf{x),Xg{x)) 

= ^c{pkAx)) {Xj{pkAx)) ,X-g{pKSx))) = {f,g}{pKAx)), 

where / = {pKc)*f: 9 = {PKc)*9- Denote by itj{c): R^^^ — )■ R^^^I^^^ the natural 
homomorphism of Poisson algebras. Then from the previous equalities we obtain 
{pkJ^Kc) {Tri(c)if),TTj(c){g)} = {f,g}, which shows that {pkJ^i(c) is a Poisson iso- 
morphism. In other words, we have the commutative diagram 

(PKc)* X y 

HMc) 

All the arrows are defined in a natural way and they are surjective Poisson algebra 
homomorphisms. ■ 
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Physical variables. It is immediate to see from the definitions that for aU x £ Mc, 

kera;(3;) n Tr,Mc C kerujc{x). (3.5) 
From now on we wih assume the following. 

Assumption 3.18. (a) keiuj(x) is a regular distribution, that is, it determines a 
regular foliation K and the natural projection px '■ M — )• M , where M = M/K , is a 
submersion. 

(b) keTUj{x) n TxMc is a distribution of constant rank. 

Theorem 3.19. The distribution kera;(j;) n TxMc is regular and has rank d^'^\x). 
Its integral manifolds are SCiMc, where S is an integral manifold ofkerco. Moreover, 
these integral manifolds give a foliation of Mc which is regular, that is, the natural 
map Pj^ : Mc — )• Mc, where Mc = Mc/Kc is a submersion. Besides, each leaf of the 

foliation Kc is foliated by leaves of Kc, which gives a naturally defined submersion 
Pj^ : Mc — )• Mc. In other words, we obtain the commutative diagram 

Mc^Mc 




where each arrow is a naturally defined submersion. 

Proof. The first assertion, about the rank of the distribution ker a;(x)nTa;Mc, is easy 
to prove. Let xq G Mc. Then there exists a uniquely determined integral manifold 
5 of the distribution kerw such that xq G S. Using that, by assumption, kera;(x) n 
TxMc is a distribution of constant dimension and that dim (kera;(x) n T^^Mc) = 
dim {TxS n TxMc) we can conclude that the intersection S n Mc coincides with the 
integral leaf of the integrable distribution of ker uCiTMc containing xq. So we obtain 
the foliation Kc of Mc- Using (3.5) we can deduce that each leaf of the foliation Kc 
is foliated by leaves of Kc. The rest of the proof follows by standard arguments. ■ 

Lemma 3.20. (a) The following diagram is commutative 
M M -~ TM 




where the arrows are defined as follows. The maps px, PKci Pk '^"■'^ Pk k 
defined in Assumption 3.8, Assumption 3.18 and Theorem 3.19. By definition, the 
map fc is the inclusion. The map fc is an embedding defined by fc{S n Mc) = S, 
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where S is a leaf of the foliation K. We will think of fc as being an inclusion. The 
vector bundle TM\Mc is the tangent bundle TM restricted to Mc- In other words, 
since fc is an inclusion, TM\Mc is identified via some isomorphism, called Cc, with 
the pullback ofTM by fc- We call Fc the natural map associated to the pullback. 

(b) The presymplectic form u on M passes to the quotient via px giving a 
uniquely defined symplectic form Co on M, satisfying p*j^u} = co. The presymplec- 
tic form LOc, which, by definition is f*uj, defines uniquely a presymplectic form 
Zjc on Mc via p^ satisfying p*~ Uc = ^c, = f*^- The energy S on M satis- 
fies d£{x)\'keTUj{x) = 0, for all x E Mc, therefore it defines uniquely a 1-form on 
TM\Mc, called {Fc o efdS G T {{TM\Mc)*^. Since £ is constant on each leaf of 

Kc, it also defines a function Ec on Mc- Since TMc C TMc via the inclusion Tfc 
we have {Fc o e)*d£\TxMc = d£c{x), for all x G Mc- 

(c) Equation of motion (3.4) on Mc passes to the quotient Mc as 

Ld{x){X{x), ) = {Fc° e)*d£{x), (3.6) 

where X{x) G T^Mc. This means that if X{x) £ TxMc is a solution of (3.4) then 
X{x) := TxPj^ X{x), where x = Pj^ {x), is a solution of (3.6). Therefore, a solution 

curve x{t) of (3.4) projects to a solution curve x{t) = {x{t)) of (3.6) on Mc. 

Equation (3.6) has unique solution X(x) for each x £ Mc- This solution also satisfies 
the equation 

ojc{x){X{x), ) = d£{x). (3.7) 

However solutions to equation (3.7) are not necessarily unique, since kera;c(x) is 
not necessarily 0. 

(d) The restriction of the energy £\Mc satisfies 

d{£\Mc){x)\kerUc{x) = 0, 

for all X G Mc, therefore there is a uniquely defined function £c on Mc such that 
PKc^c = £\Mc. The equation 

ojc{x){X{x), ) = d£c{x) (3.8) 

has unique solution X{x) for x G Mc. If X{x) is a solution of (3.6) then X{x) = 
TxPx K " solution of (3.8). Therefore, a solution curve x{t) of (3.6) projects 

to a solution curve x{t) = ^ {^{i)) of (3.8) on Mc. 

Proof, (a) The equality px^ = Pj^ ^ o Pj^ was proven in Theorem 3.19. The 
equahty T^joFcoec = fc°TM results immediately from the definitions. The equality 
Pk ° fc = fc° Px results by applying the definitions and showing that, for given 

X G Mc, Pk ° fc{x) = Sx = fc° {x), where Sx is the only leaf of K containing x. 

(b) Existence and uniqueness of a) and w is a direct consequence of the definitions 
and standard arguments on passing to quotients. For any x G Mc we know that 
there exists a solution X of equation (3.4), from which it follows immediately that 
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d£{x) \ kera;(x) = 0. The rest of the proof consists of standard arguments on passing 
to quotients. 

(c) We shall omit the proof of this item which is a direct consequence of the 
definitions and standard arguments on passing to quotients. 

(d) If X{x) is a solution of tj{x){X(x), ) = d£{x) then it is clear that it also 
satisfies uJc{x){X{x), ) = d{£\Mc){x). It follows that 

d{£\Mc){x)\kevu}c{x) = 0, 

for all X G Mc- The rest of the proof is a consequence of standard arguments on 
passing to quotients. ■ 

Remark 3.21. Recall that the locally constant function d^'^^x) on Mc is the dimen- 
sion of the distribution kera;(x) R TxMc on Mc and also the dimension of the fiber 
of the bundle S^^l U d^'^^x) is nonzero then there is no uniqueness of solution to 
equation of motion (3.4), since solution curves to that equation satisfy, by definition, 

Lo{x){±, ) = d£{x)\TxM 
(x, x) £ T^Mc. 

Passing to the quotient manifold Mc eliminates this indeterminacy and unique- 
ness of solution is recovered. This is related to the notions of physical state and 
determinism, mentioned by Dirac [1964, p. 20]. 

The reduced equations of motion on the symplectic manifold Mc are Hamilton's 
equations. The reduced equations of motion on the manifold Mc are not necessarily 
Hamilton's equations but they are ODEs. The natural fibration ^ : Mc — >• Mc 

has the following property. For rfiQ £ Mc and fho G (p^ ^ )~^{mo) C Mc, if m(t) 
and m(t) are the unique solution curves through friQ and mo respectively, then 
Pr if ™'(^) ~ 'rn{t). This means that there is a one-to-one correspondence between 

the solution curves on Mc and Mc, modulo the choice of a point on the fiber over 
the physical state tjiq. See also remark 3.30. 

3.3 Second class subbundles 

The notions of first class and second class constraints and functions and also second 
class submanifolds only depend on the final constraint submanifold Mc and the 
ambient symplectic manifold and do not depend on the primary constraint Mq = 
M or the energy £ : — )• M. Accordingly, in this section we will adopt an abstract 
setting, where only an ambient symplectic manifold and a submanifold are given. 
Among several interesting references we cite Sniatycki [1974] which has some points 
of contact with our work. 

The following definition is inspired by the ones given in section 3.2. 

Definition 3.22. Let (P, fi) be a symplectic manifold and S C P a given submani- 
fold. Then, by definition, 

r(S,p) .- {/ e ^(p) I Xf{x) G T^S, for all x e S} 

/(^'^) :={/Gi?(5'^) 1/15 = 0} 
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Elements of R^^'^^ are called first class functions. Elements of are called 

first class constraints. The submanifold S is called a first class constraint 

submanifold if for all f G ^{P) the condition f\S = ^ implies f £ that is, 

I(S,P) ^ where I{s,P) the ideal of the ring T{P) of all functions vanishing 

on S. 

Obviously, using the notation introduced before in the present section, i?(^^c,JV) — 
and = /(<=). All the properties proven for i?^^) and /^"^^ hold in general 

for and For instance, T^S is the set of all / G R^^^P\ Ev- 

ery function / G R^^'^^ satisfies df{x){Xx) = 0, for all € keT{^{x)\TxS) and 
kei{n{x)\T^S) is the set of all Xf{x), f G We shall collect in a single lemma 

3.23 the basic properties of first and second class constraints and Dirac brackets 
and their proofs in a compact way and with a uniform language that does not rely 
on coordinates. Most of these results are essentially known since the fundamental 
works of Dirac, and have been stated and proven in different ways in the literature, 
except perhaps for the notion of second class subbundle and the corresponding de- 
scription of second class submanifolds as tangent to the subbundle. The second class 
subbundles Q TP\S classify all second class submanifolds containing 5 at a 
linear level, that is, is tangent to the second class submanifold. For such a sec- 
ond class submanifold, which is a symplectic submanifold, the Dirac bracket of two 
functions F and G at points x £ S can be calculated as the canonical bracket of the 
restrictions of F and G. This has a global character. A careful study of the global 
existence of a bracket defined on sufficiently small open sets [/ C P containing S 
which coincides with the previous one on the second class submanifold will not be 
considered in this paper. However, to write global equations of motion on the final 
constraint submanifold one only needs to know the vector bundle V^, which carries 
a natural fiberwise symplectic form. 

All these are fundamental properties of second class constraints and constraint 
submanifolds, and Theorem 3.24 collects the essential part of them; we suggest to 
take a look at it before reading Lemma 3.23. 

Lemma 3.23. Let (-P, ^2) be a symplectic manifold of dimension 2n and S Q P a 
given submanifold of codimension r. Let to be the pullback ofVt to S and assume that 
kercL'(x) has constant dimension. Assume that S is defined regularly by equations 
(pi = 0, . . . , (j)r = on a neighborhood U ^ S and assume that we can choose a 
subset {xi, . . . , X2s} ^ {01, • • • 1 (pr} such that det({xj, Xj}{x)) /o'^ '^^^ x £ U, 
where we assume that 2s = icank{{(pi,(j)j}(x)), for all x G U. We shall often denote 
c^j{x) = {XiiXj}{x) O'f^d c-l{x) the inverse of c^-{x). Moreover, we will assume that 
the following stronger condition holds, for simplicity. Equations (pi = Bi, . . . ,(j)r = 
Br and xi = Gi, . . . , X2s = G2s define submanifolds ofU regularly, for small enough 
Bi, . . . , Br and Ci , . . . , C2s ■ Then 

(a) 2s = r — dimkercj = 2n — dim 5 — dimkerw. There are G L'^^'^\ k = 
1, . . . ,r — 2s, which in particular implies {ipk,'4'l}{x) = 0, {ipk, Xi}{x) = 0, for k, I = 
1, . . . ,r — 2s, i = 1, . . . ,2s, such that (i^i(x), . . . , dipr~2six), dxi{x), . . . , dx2s{x), are 
linearly independent for all x £ S. Moreover, X^-^(x), . . . ,X^^_2^(x) form a basis of 
kerc<j(x), for all x G S and dipi{x), . . . ,dipr-2s{x),dxi{x), . . . ,dx2s{x) form a basis 
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ofiT^sy. 

(b) The vector subbundle C TP\S with base S and fiber 

= span {X^, {x),..., X^,^ {x)) C T^P, 

satisfies 

V} n T^s = {0} 
(v;^)^n(kerw(x))^ = r,5, 

xeS. 

(c) There is a neighborhood U of S such that the equations Xi = 
on U define a symplectic submanifold such that S C. and 

T^S^ = {V^f 
T^S^ © = T^P, 

for x G S^, where we have extended the definition of Vx for x & using the 
expression 

= span {X^^ (x),..., X^,^ {x)) C T^P, 

forx G S^. The submanifold has the property I(^s,S^) ^ that is, S is a first 

class constraint submanifold of S^, defined regularly by ipi\S^ = 0, i = 1, . . . , r — 2s, 
and tpilS^ £ ^ i = 1, . . . ,r — 2s. Moreover, it has the only possible dimension, 

which is dim 5^ = dim/S + dimkerw = 2n — 2s, for symplectic submanifolds having 
that property. It is also a minimal object in the set of all symplectic submanifolds 
Pi ^ P, ordered by inclusion, satisfying S Pi. 

(d) Let V be any vector subbundle ofTP\S such that 

V®kevoj = {TS)^, (3.12) 

or equivalently, 

© (kerw)^ = {TSy. (3.13) 

Then dimV^ = 2s, for x ^ S. Let be a submanifold of P such that T^S^ = V^, 
for each x & S. Then S is a submanifold of . Such a submanifold always 
exists. Moreover, for such a submanifold there is an open set U C. P containing S 
such that nU is a symplectic submanifold of P. 

Let X £ S and let x'^ = 0, . . . , X2s = be equations defining H U' for some 
open neighborhood U' P and satisfying that dx'iix), . . . ,dx2s{x) are linearly in- 
dependent for X G n U' . Then, dxiix), . . . , dx2s{^) i dipi{x) , • • • ' dtpr-2six) are 
linearly independent and det{{xi, x'j}{x)) 7^ 0? for x G fl U' . All the properties 
established in (a), (h), (c) for xi, ■ ■ ■ , X2s on S hold in an entirely similar way for 
x'l , . • . , X2s ' on S nU' . In particular, H U' = . 

(e) Let uj^ be the pullback of to and {,}x the corresponding bracket. For 
given F,G e T{P) define F^^ := F — XiCx{Xj->P} ^^rid also 

:= {F,G} - {F,Xi}4{X3,G}. 



(3.9) 
(3.10) 
(3.11) 

, • • ■ , X2s = 0, 
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which is the famous bracket introduced by Dirac, called Dirac bracket, and it is 
defined for x in the neighborhood U where cfj (x) has an inverse (x) . Then, for 
any x G S^, 

{F^,Xi}(x) = 

for i = 1, ... ,2s, and also 

{F^,G^}{x) = {F,G}^^^{x)=uHx) (X,.|5.(.t),Xg|sx(x)) = {F\S^,G\S^^{x). 

If we denote the Hamiltonian vector field associated to the function F G 

T{P), with respect to the Dirac bracket {,}{x) then the previous equalities are equiv- 
alent to 

The Jacobi identity is satisfied for the Dirac bracket {F, on S^, that is, 

{{F,G}(^),i7}(^)(x) + {{i?,F}(^),G}(^)(x) + {{G,iJ}(^),F}(^)(x) = 0, 
for X £ S^. 

(f) Let U be an open neighborhood of S such that cfj{x) is invertible for x £ U. 

For each C = (Ci, . . . , C2s) £ 1^^'^ let xf = Xi~ Ci and define S^^ by the equations 

Xp(x) = 0, i = l,...,2s, X £ U. For any C in a sufficiently small neighborhood 

c ^ 
of 0, is a nonempty symplectic submanifold of P. Define the matrix cfj (x) = 

{xf\x9'}{x) , and also c-'dx) as being its inverse, x £U. Then, the equalities 



cfj{x) = {Xi,Xj}{x) = {x?,Xj}{x) = cf/(x), (3.14) 

and also, 

{F,G}(^)(x) = {F,G}(^c)(x) (3.15) 

are satisfied for all x £ U. All the definitions and properties proved in (e) for the 
case C = hold in general for any C in a neighborhood of small enough to ensure 
that is nonempty. In particular, the equalities 

{F^c,G^c}{x) = {F,G}^^c){x) =uj^''{x) (Xj,|^^c(x),X^|^^c(x)) 

= {F\S^'',G\S^''}^c{x) (3.16) 

and 

hold for X £ , and any C in such a neighborhood. The Dirac bracket {F, G}(^) 
satisfies the Jacobi identity for F,G £ J^{U) and the symplectic submanifolds 
are the symplectic leaves of the Poisson manifold {F{U),{, }(x))- shrinking, if 
necessary, the open set U and for C in a sufficiently small neighborhood ofO£ M?^, 
the equations ipk\S''^ =0; k = 1, . . . ,r — 2s, define regularly a first class constraint 

submanifold C S^^ C U, and the functions il^klS^^ e Ri^^'^'''') C J^{S^^) are 
first class constraints, that is, ijjk\S^'^ £ l(^^'^^ ) k = \, . . . ,r — 2s. We have that 
dimS'^ = dimS''^ + dimker w'^, where uj^ is the pullback of 17 to 5*^. One has 
dimS''-^ = dim 5* and dim 5^ = dimS"^, therefore dimkerw = dimkerw*^. 



3.3 Second class subbundles 23 



Proof, (a) Let x £ S. We are going to use lemmas and corollaries 3.1-3.7 with 
E := T^P; V := TxS; 7j := d4)i{x), z = 1, . . . , r; O := oj := oj[x); /3 = 0. 

Elements ip = X'(j)i, A* G ^(P) such that V' G which implies G TxS 

for x G 5, must satisfy {'tp,(pj}{x) = 0, or, equivalently, A*(x)d(/)i(x) = 
0, for j = l,...,r, x G 5". Using Lemma 3.6 we see that X^(x) G keT:u){x). 
Since one can choose r — 2s linearly independent solutions, say Aj = (Xj, . . . , A[), 
i = l,...,r — 2s, we obtain elements ipi G namely, ipi = \l4>j^ such that 

{dil)i{x), . . . ^dil)r-2s{x)) are linearly independent, or, equivalently, taking into ac- 
count Corollary 3.7, that {X.ip^{x), . . . , X^^_^^{x)) is a basis of kercj(x) for x £ S. 
If d^i{x), . . . , dipr-2s{x), dxi{x), . . . , (ix2s(2;), were not linearly independent, then 
there would be a linear combination, say x = CL^Xi^ with at least one nonzero coeffi- 
cient, and some x £ S, such that dx{x) = n'^dipkix) for some fi^, k = 1, . . . ,r — 2s. 
But then, for any j = l,...,2s, {x,Xj}{x) = dx{x)X^^{x) = fi''dipk{x)Xy^^{x) = 
—lJL^dxj{x)X^^{x) = 0, which contradicts the fact that det {{Xi^Xj}{x)) 7^ 0. Using 
this and the fact that ^/^j = 0, = 0, i = r — 2s, j = 1, . . . , 2s define S regularly, we 
can conclude that dijji{x), . . . , dipr-2s{x), dxiix), . . . , dx2six) is a basis of (TxS)°. 

(b) If X'Xy.^{x) G TxS then dxj{x)X'X^^{x) = 0, j = 1, . . . , 2s, which implies 
^^{Xji Xi} = 0, j = 1, . . . , 2s, then A* = 0, i = 1, . . . , 2s, which proves (3.9). To 
prove (3.10) we apply the operator to both sides and obtain the equivalent equality 
span((ixi(3;), • • . ,dx2s{x)) ® span{d7pi{x) , . . . , di/jr-2s{x)) = (TxS)°, which we know 
is true, as proven in (a). To prove (3.11) we apply the orthogonal operator ^ to 
both sides of (3.10). 

(c) Since dxi{x), . . . ,dx2s{x) are linearly independent for x G 5 they are also 
linearly independent for x in a certain neighborhood U of S. Then xiix) = 
0, . . . ,X2s{x) = define a submanifold S^ of U containing S. To see that it is 
a symplectic submanifold choose x £ S^ and apply Corollary 3.7 with E := TxP; 
V := TxS^; j3 := 0; 7j := dxi, i = I, ■■.,2s; oj := VL{x)\TxS'^^ We can con- 
clude that dim(ker r2(x)|T^5'^) = 0. Now let us prove that TxS^ = {Vx)^, namely, 
TxS^ = span(dxi(x), . . ■,dx2s{x)f = {{V^ff = {V^f^ From this, using that S^ 
is symplectic one obtains TxS^ ® Vx = TxP^ To prove that 5" C 5-^ is a first class 
constraint submanifold defined by first class constraints 'ipi\S^, i = 1, ■ ■ ■ ,r — 2s, on 
S^, we observe first that it is immediate that ■i/'il'S'^ = 0, i = 1, . . . , r — 2s, define S 
regularly. It remains to show that X^.\gx{x) G TxS, i = 1, ■ ■ ■ ,r — 2s, x £ S , where 
X^.^gx is the Hamiltonian vector field associated to the function ipi\S^ with respect 
to the symplectic form w-^. This is equivalent to showing that 

x^^\sAx) {i^j\s^) = o, 

for X G 5 or, equivalently, 

u^{x) (^X^^|5x(x),X^^.|5x(2:)) = 0, 

for X £ S. We know that uj^ is the puUback of to S^ and ipilS^ is the pullback of 
il^i, via the inclusion S^ U, then we have 

uj^{x) (^X^^|5x(a;),X^,|5x(x)) = n{x) {X^^{x),X^^{x)) = 0, 



3.3 Second class subbundles 24 



for X £ S, since ipi are first class constraints, i = 1, . . . ,r — 2s. Finally, using 
the definitions we can easily see that dim = 2n — 2s and that dim S = 2n — r 
and from (a) we know that dimkerw = r — 2s. We can conclude that 2n — 2s = 
dim S + dim ker w. 

(d) We know that for x £ S, dim T^^S" = 2n — r, and dimkera;(x) = r — 2s; then 
using (3.12) we obtain dim 14 = 2s. Also from (3.12) we immediately deduce by 
applying ^ to both sides, 

nikerojf = TS, 

in particular TS C V^. Let g he a given Riemannian metric on P and let Wx be 
the (^-orthogonal complement of TxS in V^, in particular, Wx ® TxS = V^, for each 
X G S. Define 

= {exp{twx) I Wx e Wx,gix){wx,Wx) = l,\t\ < t{x),x £ S}. 

By choosing t{x) appropriately one can ensure that is a submanifold and, more- 
over, it is easy to see from the definition of that TxS^ = Wx ® TxS = V^, for 
each X £ S. We leave for later the proof that n C/ is a symplectic submanifold 
of P, for an appropriate choice of the open set U, which amounts to choosing t{x) 
appropriately. 

Assume that dxi{x) , . . . , dx2six) are linearly independent for x E fl U'. 
Since {dXi{x),V^) = {dXi{x),TxS'^) = for x S and i = l,...,2s, we can 
deduce that dx'i{x) € {V^Y , that is, dx'iix) £ Vx- Then using (3.13), we see 
that dx'iix)., . . . , dx'2s{x) 1 dx'iix) , . . . , dipr-2six) are linearly independent and span 
y^®(kerw)^ = {TSy. If det({x-, Xj}(^)) = for some x e 5 then A*{x-, X^Ka^) = 0, 
where at least some A* 7^ 0, i = l,...,2s. Let X^x'i = then {ip,x'j}ix) = 0, 
j = l,...,2s. On the other hand, since (p\S = 0, then {ip,ipj}{x) = 0, j = 
l,...,r — 2s. We can conclude that if G /(■^''^ ) and then Xip{x) G kertj(x), in 
particular, X^{x) = fi^X^_^{x), which implies X^dx'iix) = dip{x) = fiMipjix), contra- 
dicting the linear independence of dx'i{x), . . . , dx'2s{x), dilji{x), . . . , diljr-2s{x). 

It follows from which precedes that by replacing Xi by x'j) ^ = 1> • • • > 2s and S 
by 5 n [/' all the properties stated in (a), (b) and (c) are satisfied. In particular, 
S nU' = and S CiU' is symplectic. It is now clear that by covering S with open 
subsets like the U' we can define U as being the union of all such open subsets and 
one obtains that S CiU is a symplectic submanifold. 

(e) Let X G S^. Then, since = F — Xi^xiXji P}: we obtain 

{Fx,Xk}{x) = {F,Xk}ix)-{xi,Xk}ix)d^ix){Xj,F}ix) = {F,Xk}ix)+{Xk, F}{x) = 0. 
Using this we obtain 

{F^, G^}{x) = {F^, G - Xkc'^ixi, G}}{x) = {F^, G}{x) 

= {F,G}(x) - {x^,G}c^{xJ,F} = {F,G}(^)(x). 

For any F £ F{P), x £ and A; = 1, . . . , 2s, we have Xp^ ix)xk = {Xk, Fx}i^) ~ 
0, so Xp^{x) £ TxS^. Therefore, for any Yx £ TxS^, 

w^(x) {XF^ix),Yx) = nix) {XfJx),Yx) = dF^ix)Yx = diF^\S^) ix)Yx 

= diF\S^) ix)Yx = io^x) {Xp^s^ix),Yx) , 
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which shows that Xp^{x) = Xp^gxix), where both Hamiltonian vector fields are 
calculated with the symplectic form uj^. Using this, for any G G J~{P) and any 
X £ S^, one obtains 

{G^,Fj(x) = Xf^{x)G^ = Xf^{x)G = Xf\sAx)G\S^ = {G\S^ , F\S^ ^{x) . 

The equality X(^^^^p{x) = Xf^{x) = Xp\sx{x) is an immediate consequence of the 
previous ones. The Jacobi identity for the bracket { , }(-^) follows using the previous 
formulas, namely, for x G S-*^, one obtains 

= {{F,G}^^p^,H\S^}^{x) = {{F\S^,G\S^^,H\S^Ux), 

where the bracket in the last term is the canonical bracket on the symplectic manifold 
S^, for which the Jacobi identity is well known to be satisfied. 

(f) The equalities (3.14) and (3.15) are proven in a straightforward way. The 
equations (3.16) and (3.17) follow easily using a technique similar to the one used 
in (e). Using all this, the proof of the Jacobi identity for the bracket { , j^^^,) on [/ 

goes as follows. Let x £ U and let C be such that x £ . For F,G,H G F{U) 
using (e) we know that the Jacobi identity holds for { , }(j(.c-) on . But then, 
according to (3.15) it also holds for {, }(^) for all x £ . Now we will prove 
that S^^' are the symplectic leaves. Since they are defined by equations xp = 0; 
i = 1, . . . ,2s onU we need to prove that {F, xf}x(^) = 0, x G S*^, for all F G F'{U), 
i = 1, ... ,2s. Using (3.15) and (3.16) we see that {F, xF}(x)(^) — xf}{x'^){x) = 

,'XX\^^ }x'^(^) ~ ^- finish the proof, observe first that, since Xi = 0; 
V'i = 0, i = 1, . . . , 2s, j = 1, . . . , r — 2s define regularly the submanifold 5 C {/, by 
shrinking \J if necessary and for all G sufficiently small, we have that = 0, ^/^j = 0, 
i = \,. . . , 2s, j = 1, . . . , r — 2s define regularly a submanifold S'~'' C U and therefore 
ipj\S^ = 0, j = 1, . . . , r — 2s, define regularly 5*^ as a submanifold of . To prove 
that it is a first class constraint submanifold and that i^jlS^ , j = 1, . . . , r — 2s, 
are first class constraints, that is ipj\S^^ G I^^^'^^ \ j = 1, . . . ,r — 2s, we proceed 
in a similar fashion as we did in (c), replacing X by X • The fact that i/^jlS^ , 
j = 1, . . . ,r — 2s, are first class constraints defining S'^ implies that dimS'^^'^ = 
dim 5*-^ + dimkerw*-". From the definitions one can deduce that dim 5*-" = dim 5" 
and dim S"-^ = dim 5-^, therefore dim kerw = dim ker a; ■ 

The following theorem summarizes the essential part of the previous lemma. 

Theorem 3.24. Let {P, 17) be a symplectic manifold, S" C P and letuj he the puUback 
of Vt to S. Assume that kerw has constant rank. Let V he a vector subbundle of 
TP\S such that F kerw = (TS)^. Then there is a symplectic submanifold 
containing S of dimension dim 5" + dimkerw such that the condition T^S^ = , 
for all X £ S holds. The vector bundle is called the second class subbundle 
tangent to the second class submanifold . The submanifold is locally defined 
by certain equations Xi = 0; ■ • • > X2s = as in Lemma 3.23. Under the hypotheses 
of that lemma, namely, that 2s is constant on a neighborhood U of S, will be 
one of the symplectic leaves of the Dirac bracket { , }(x)- 
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Remark 3.25. (a) Under our strong regularity conditions the symplectic leaves of the 
Dirac bracket give a (local) regular foliation of a neighborhood of the final constraint 
submanifold S. This implies by the Weinstein splitting theorem (Weinstein [1983]) 
that there are local charts where the Dirac bracket is constant. 

(b) The tangent second class subbundle in a sense (modulo tangency) clas- 
sifies all the possible second class constraint submanifolds containing a given sub- 
manifold C P. It carries enough information to write the Dirac brackets along the 
final constraint submanifold S and therefore also equations of motion, as we show 
in section 3.4. 

3.4 Equations of motion 

We are going to describe equations of motion in the abstract setting of section 3.3, 
that is, a symplectic manifold {P, and a submanifold S C P, defined regularly 
by equations (pi = 0, i = 1, . . . , a. For this purpose, we need to introduce in this ab- 
stract setting, by definition, the notions of primary constraints, primary constraint 
submanifold and the energy. 

The primary constraint submanifold is a given submanifold S" C P con- 
taining S, and in this context, S will be called the final constraint. We will assume 
without loss of generality that S' is defined regularly by the equations (pi = 0, 
i = 1, . . . ,a' , with a' < a, where each (pi, i = 1, . . . ,a' will be called a primary 
constraint while each (pi, i = a' + 1, . . . ,a will be called a secondary constraint, 
for obvious reasons. In this abstract setting the energy is by definition a given 
function £ G J"(P). 

The equations of motion can be written in the Gotay-Nester form, 



where {x,x) G T^S, for all dx G T^S' . We require that for each x G S, the space of 
solutions {x,x) of (3.18) is nonempty and forms an affine bundle. 

Now we will transform this equation into an equivalent Poisson equation using 
the Dirac bracket. 

The condition {£,ip}{x) = 0, for all x S 5* and all first class constraints ip will 
appear later as a necessary condition for existence of solutions for any given initial 
condition in S, so we will assume it from now on. 

The total energy is defined by 



where the functions A* G C°°(P), i = 1, . . . , a' must satisfy, by definition, {£t, 'Pj}ix) = 
0,j = l,...,a,x£Sor, equivalently, 



^l{x){x,5x) = dS{x){5x) 



(3.18) 



£t = £ + ><(pi, i = l,...,a' 



(3.19) 



{£, (t)j}{x) + X'{(j)i, ^j}{x) = 0,xeS,i 
sum over i = 1, . . . ,a' . 



= l,...,a,j = l,...,a. 
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Assumption 3.26. We assume that the solutions (A"^, . . . , A"') form a nonempty 
affine bundle A — t- S*. Equivalently, if we call d\ the dimension of the fibers, s' := 
rank({0i, = a' — d\ is constant on S. 

By taking a sufficiently small neighborhood U, we can assume that A is a trivial 
bundle over S CiU, which we extend arbitrarily as a trivial affine bundle on U. 

For each section of A one has a Hamiltonian £t{x) = £{x) + A*(x)(/»j(x), x e U, 
and the equation of motion on S CiU, 

Xs^ = {dSrf. 

Equations of motion in terms of the Dirac bracket. Local equations of 

motion in terms of a Dirac bracket become simpler when the second class constraints 
are adapted to the primary constraint, as we will explain next. This has been 
considered by Dirac. The case of general second class constraints is important for 
the global case, and will not be studied here. 

Assume that the primary and final constraints S' D S are defined by 

0. = 0, i = 1, . . . , a', defines S' (3.20) 
01 = 0, i = 1, . . . , a', a' + 1, . . . , a, defines S, (3.21) 

regularly on a certain neighborhood U of 5". For choosing first class and second 
class constraints one can proceed as follows. 

Assumption 3.27. The rank of the matrix $(x) = {{(pi,(pj}{x)), i,j = l,...,a, 
which we will call 2s, is constant on an open neighborhood U of S' . 

By shrinking U if necessary, one can choose {xi, . . . , X2s} ^ {^ij ■ ■ ■ , 4'a} in 
such a way that ({%?, Xj}(ic)) is a regular submatrix. Without loss of generality, 
assume that these functions are the last 2s, that is, Xi = 4'a-2s+i, i = 1, • . • , 2s. For 
each (pi, i = I, . . . ,a — 2s, there are af{x) such that ijji{x) := 4>i{x) + a^(a;)x/;(a;), 
A: = 1, . . . , 2s, defined on U satisfy 

{ipi, (t)i}{x) = {(pi, (j)i}{x) + (x){xfe, Mix) = 0, 

for all x e S, and alH = 1, . . . , a. Therefore, these ipi are first class constraints and 
{Xij • ■ ■ 5 X2s} are second class constraints, and for each x e U the 

dipiix), dipa-2s{x),dxi{x), dx2s{x), 

are linearly independent, x £ U. We will sometimes call (ipi, . . . tpa-2sj Xij • • ■ X2s) a 
complete set of first class and second class constraints for S. 

Now we will show how to choose a complete set of first class and second class 
constraints adapted to the primary constraints. This, as a definition, means 

choosing 

Iplix), 1pa-2s{x), Xl{x), X2s{x), 
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as before in such a way that, for s' equal to the rank of the submatrix = 
({<^ii^j}W)> i = l,-..,a',j = l,...,a, 

V'l,...,V'o'-s',Xl,---,Xs' 

are primary constraints. 

To achieve this, let {x'l, • • • , x's'} ^ {4*1^ • • • > 4'a'} such that the rows {x^, (j)j}{x), 
k = 1, . . . ,s' , are linearly independent. Without loss of generality, assume that 
these functions are the last s' , that is, x'i = <Pa'-s'+i, i = 1, ■ ■ ■ , s' . Then, for each 
i = 1, . . . ,a' — s', there are (x), k = 1, . . . ,s' such that 

4{x) := cl)i{x) + al{x)x'k{x) 

defined on U satisfy 

{4, Mix) = Mix) + a'^ix){x'k, Mix) = 0, 

for all i = 1, . . . , a' — s', x G S, and all I = 1, . . . , a. Therefore, these ipi form 
a set of a' — s' first class constraints. Note that (-01, . . • V'a'-s'' x'l' • • • > ^s') ^ 
complete set of primary constraints, in the sense that their differentials are linearly 
independent on U, and every primary first class constraint ip can be written as 
ipix) = biix)'tp[ix) H h ba'-s'ix)tp'a,_g,ix). 

In order to obtain a complete set of first class and second class constraints 
we first choose {x's'+ii • • • > x'L} ^ {4>a'+i, ■ ■ ■ , (f>a} such that the rows {x'k'-'Pj}{^)^ 
{Xfc") 'Pj}{x) k' = 1, . . . , s' , k" = s' + 1, . . . ,2s are linearly independent. This can 
be done because the matrix ({^j, i,j = l,...,a, has rank 2s. As before, 

assume that Xi = ^a-2s+i, i = s' + 1, . . . ,2s. 

Then, for each i = a' — s' + 1, . . . , a — 2s there are (x), (x), /c' = 1, . . . , s', 
k" = s' + l,...,2s, such that ^j-ix) := (ps'+iix) + af (x)Xfe/(x) + af (x)Xfc//(a;), 
k' = 1, . . . , s' , k" = s' + 1, . . . , 2s, satisfies 

{ij'/,Mix) = {(l)s'+i,Mi^') + <^i\x){Xk'Al}{x) + afix){xk»,(t>i}ix) = 0, 

for all X & S, and alH = 1, . . . , a. 

It can be shown using the definitions that 

dx'lix), dx's'ix),dxs'+iix), dx2s{x)) 

are linearly independent for x £ S. Note that no nontrivial linear combination of 
the V'f is a primary first class constraint. 
Therefore, 

(V'l,...,^a-2s,Xl,---,X2s) : = 

ilp'l,...,tpa'-s'^'^'i'-s'+l^--- 5 V'a-2s5 X'li • • • 5 Xs'j Xs'+l) ■ ■ • ) X2s) 

form a complete set of first class and second class constraints adapted to the primary 
constraints. 
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We will now take advantage of these adapted constraints to write the equations 
of motion in a simpler way. The total energy defined on U is 

fT = f + a'v: + m'^x;-, 

z = 1, . . . , a' — s', J = 1, . . . , s'. As usual, we impose on £t the condition that it is 
first class, which means that {iS, il:i\{x) = 0, z = 1, . . . , a — 2s and 

i = l,...,2s, X G S. From this we obtain well-defined n'^ix), x G S, which we 
extend arbitrarily for x G U. Then the equations of motion are written as 

X £ S, where A'*, i = 1, . . . , a' — s' , are completely arbitrary real numbers. We can 
also write this equation in the Poisson form 

F = {F,£ + fi'^x'j} + ^''{F,i^l}. 

Now we will see a simpler expression for the equations using Dirac brackets 
rather than the canonical bracket. Since £t is first class, then X(;j,{x) = £-^,(2;) 
for X £ S (see Lemma 3.23(e)). Then the equations of motion can be written as 

F = {F,£T}(^xy (3.22) 

This equation, obtained in an abstract setting, is similar to the one in page 42 of 
Dirac [1964] for the case P = T*Q. 
We can rewrite the total energy as 

£t = £ + A'>^ + /i'^x;- =£ + A'^(0. + a^x'k) + l^''Xj, 
where i = 1, . . . , a' — s' and j = 1, . . . ,s' , and therefore (3.22) can be written as 

F = {F,£}^^) + X'^{F,cP,}^^y 

or equivalently, 

xeS. 

We will denote 

£at = £ + (3.23) 

i = 1, . . . , a' — s' , and call it the abridged total energy. It has the property that it 
provides the same dynamical information as the total energy in terms of the Dirac 
bracket. Moreover, for each x £ S the map ~^ — )• given by (A'"^, . . . , A'*^ ~^ ) 1— )• 
■X(x),£at(^) is an isomorphism. In particular, for any solution curve x{t) of on 
S, we have 

{x{t),x{t)) = X(^^ysAA^{t)) 

for a (uniquely determined) choice of A'*(t), i = 1, . . . , a' — s' . 
We have proven the following theorem. 
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Theorem 3.28. Let the primary S' and final constraint S be submanifolds de- 
fined regularly by equations (3.20), (3.21) on a neighborhood U P of S' , where 
(P, ri) is a symplectic manifold. Let 2s be the rank of the matrix $ = {{(j)i,(j)j}{x)), 
i,j = l,...,a, which, by Assumption 3.27, is constant on U. Let s' be the rank of the 
matrix ^' = {{(pi, (j)j}{x)), i = 1, . . . ,a' ,j = 1, . . . ,a, which, by Assumption 3.26 is 
constant on S. Choose s' linearly independent rows ({x^i (t^j}{x)), i = 1, . . ■ , s' , j = 
1, . . . ,a of^', say w.l.o.g. x'i = 4'a'-s'+i> i = 1, ■ ■ ■ , s' . Then for each i = 1, . . . , a'—s' 
one can find coefficients a'l{x) such that i(i[{x) := (t)i{x) + a^(x)x'^(x), k = 1, . . . , s' , 
defined on U constitute a maximal set of linearly independent primary first class con- 
straints. Then {ip[, . . . , ip'^/_g/,Xii ■ ■ ■ iX's') ^ complete set of primary constraints 
and one can extend it to a complete set of first class and second class constraints 

{'ljjl,...,iJa-2s,Xl,---,X2s) ■■ = 

I 1 1 1 1 I II I II I I II II \ 

[Vli • • • ) Va'-s'lVa'-s'+li • • • ' Va-2si Xli • • • i Xs'i • • • ) X2s)- 

Let £ : U ^ be a given energy satisfying {£,'ipi}{x) = 0, xGS,i = l,...,a — 2s. 
Consider the following abridged total energy 

£at = £ + A'Vi, 

i = 1, . . . ,a' — s' , where \[ are arbitrary real parameters. The corresponding Hamil- 
tonian vector field with respect to the Dirac bracket is 

X(^),^^(x) = + A'*X(^),<^^(x), (3.24) 

X £ U. For x G S, this vector field is tangent to S and gives the equations of motion. 
Equivalently, we can write the equations of motion as 

F = {F,£:}(^) + A'HF,</'.}W. (3.25) 

Since X(^) is tangent to S, the evolution of a function f on S is given by (3.25) 
for any F such that f = F\S. 

Remark 3.29. Equation (3.25) represents an elaboration on the equation (3.22). It 
has interesting advantages which may be useful in practice, because it only requires 
calculating (xi, • • . , x's'j x's'+i' • • • ' ^2s)' which is a set of second class constraints 
adapted to the primary constraint. It will be generalized to extend its applicability 
to examples where the constraints are foliated, as we will see in section 7. 

Remark 3.30. The geometric meaning of s' is given by the fact that a' — s' is the 
number of independent primary first class constraints ip'^, i = 1, . . . ,a' — s'. This 
number coincides with d^'^'^{x) and with the dimension of the leaves of the foliation 
Kc of Mc (see Theorem 3.19). Also, the Hamiltonian vector fields X,^', . . . ,X^'^ ^ 

are tangent to the leaves of Kc. Using the formula = X_^f gy it is easy to 

prove that the iterated Poisson brackets of the tp'- give functions whose Hamiltonian 
vector fields are tangent to the leaves of Kc- Each one of these functions, let us 
call it /, satisfies that Xf{x) = X^{x)X^i{x) for all x G Mc, but not necessarily / 
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is a linear combination of the ip^, on a neighborhood of Mc, so it might not be a 
primary first-class constraint. One can prove that such a function is zero on Mc, 
but since df{x) = X^[x)dil)[{x) for all x G Mc, it can not be included in a set of 
constraint functions defining Mc regularly and containing the ^/j^, in other words, 
it is not an independent secondary first class constraint. There are examples in 
the literature where the Poisson bracket of two primary first class constraints is 
an independent secondary first class constraint, for instance, that is the case in 
Henneaux and Teitelboim [1992, p. 2]; for that example Assumption 3.18 (a) does 
not hold (but 3.18 (b) does). 

According to Dirac [1964], pages 23-24, primary first class constraints and their 
iterated brackets represent transformations that do not change the physical state of 
the system, which in this case is implied directly from the requirement of determin- 
ism. We have seen that under our regularity conditions, those iterated brackets do 
not really add any new such transformations besides the ones given by the primary 
first class constraints. 

In addition, note that a— 2s is also the dimension of the kernel of the presymplec- 
tic form on S (see Theorem 3.23 (a)), which is the number of first class constraints, 
and also the dimension of the leaves of (see Assumption 3.8). Since the Hamil- 
tonian vector fields corresponding to the first class constraints generate the 
integrable distribution associated to Kc, then the transformations that preserve the 
leaves of Kc are those which preserve the physical state. So, points of Mc = Mc/Kc 
would represent exactly the physical states. 

Dirac [1964, p. 25] introduces the notion of extended Hamiltonian. In our con- 
text, we should define the notion of extended energy, 

£e = £t + X"'i^i, (3.26) 

i = 1, . . . , a — 2s — {a' — s') where the ■0"* represent the secondary first class con- 
straints. Both £e and £t give the same dynamics on the quotient manifold Mc- 
If {Mc,u}c) is symplectic, there are only second class constraints among the (pi, 
i = 1, . . . ,a, and conversely. In this case, £e = £t- 

4 Dirac structures 

The rest of the paper is devoted to generalizing the Dirac and Gotay-Nester theories 
of constraints. In this section we review a few basic facts related to Dirac structures 
and Dirac dynamical systems following Courant [1990] and Yoshimura and Marsden 
[2006a,b], which are of direct interest for the present paper. 

We introduce the flat and the orthogonal ^ operators with respect to a given 
Dirac structure, which is important to describe the algorithm, in the next section. In 
fact, we obtain a technique which imitates the Gotay-Nester technique, by replacing 
uj by D. 
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Dirac structures on vector spaces. Let V be an n-dimensional vector space 
and let V* be its dual space. Define the symmetric pairing {{■ , •)) onV(BV* by 

{{{vi,ai), (^2,02))) = {ai,V2) + (02,^1) , 

for {vi,ai), (w2j «2) G V ®V* where (• , •) is the natural pairing between V* and V. 
A linear Dirac structure on V, (sometimes called simply a Dirac structure on V), 
is a subspace D V (B V* such that D = , where is the orthogonal of D 
relative to the pairing ((•,•)). One can easily check the following result. 

Lemma 4.1. A vector subspace D V ® V* is a Dirac structure on V if and 
only if it is maximally isotropic with respect to the symmetric pairing. A further 
equivalent condition is given by diuiD = n and (((ui, ai), (v2, 02))) = for all 
{vi,ai), {V2,a2) e D. 

One of the main examples is the Dirac structure naturally defined on a presym- 
plectic vector space (V, ui) by 

= {{v,a) eV®V* \ a = J{v)}. 

Here the definition of J is the standard one, namely, uj^ : V ^ V* is defined by 
ljJ^{v){w) = u{v,w), for all v,w € V. We also recall the definition of orthogonality 
on V associated to a given presymplectic form u. For a given subset W V we 
define the w-orthogonal complement by 

W"^ = {v€V\ Lo{v, w) = for aU w G W}. 

We recall some standard facts in the following lemma, omitting the proof. 

Lemma 4.2. Let lo be a presymplectic form on a vector space V and letW^V 
be any vector subspace. Then := uj^{W) = [W^Y , where the right hand side 
denotes the annihilator ofW^. 

The following proposition is a direct consequence of propositions 1.1.4 and 1.1.5 
in Courant [1990]. 

Proposition 4.3. Given a Dirac structure D C V®V* define the subspace Eo ^ V 
to be the projection of D on V. Also, define the 2-form ud on Ed by uj£)(v,w) = 
a{w), where v ® a £ D. (One checks that this definition of ujd is independent of 
the choice of a). Then, ud is a skew form on Ed. Conversely, given a vector space 
V , a subspace E <^ V and a skew form u on E, one sees that D^^ = {v (B a \ v £ 
E, a{w) = uj{v, w) for all w € E} is the unique Dirac structure D on V such that 
Ed = E and ujd = 

D is the Dirac structure associated to a presymplectic form ui on V , as explained 
before, if and only if Ed = V and to = ujd- 
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The operators D and . There is a natural extension, which is important in 
the present paper, of the previous definition of uJ' for the case of a general Dirac 
structure D C V ®V* . For a given Dirac structure D ®V* and given X G E£, 
define the set sometimes denoted for short, by 

D\X) = {a eV* \ {X,a) := X ® a e D}. 

Note that D is the Dirac structure associated to a presymplectic form on V, that 
is, D = Di^, if and only if = V and for each X ^ V, the set D^{X) has a single 
element, more precisely, D\X) = {J{X)}. In this sense, D\X) generalizes J{X). 
For a given subset W CEd define D\W), also denoted W\ by 

d\w)= \J d\x). 

If is a subspace, then D^{W) is a subspace of V* . 
It is straightforward to check that for all X, y G Eq, 

{ud{X,Y)} = d\x){Y) = -D\Y){X) = -{ud{Y,X)}. 

The notion of orthogonal complement with respect to uj can be generalized as 
follows. For any subset W C E^ define by 

= {X eV\ D\Y){X) = {0}, for ah Y G W}. 

Clearly, 

HEd = {X e Ed\ d\y){X) = {0}, for all Y G W}. 
It is easy to check that for any subspace W C Ed one has 

HEd = W'^^. 

We recall also that, since ujd- Eq x Ed — t- M is a presymplectic form on Eq, one 
has, according to Lemma 4.2 and with a self-explanatory notation. 

The following proposition generalizes Lemma 4.2 and is one of the ingredients of the 
constraint algorithm described in section 5. 

Proposition 4.4. Let D (^V ®V* he a given Dirac structure and let W C Ed he 
a given suhspace. Then D\W) = = 

Proof. We first show that C {W^)°. Let a G W\ say a G D\X) for some 
X eW. Since D\X){Y) = {0} for all Y e by definition of W^, it follows that 
a G {W^y. To prove the converse inclusion, we first observe that it is equivalent 
to prove that (1^^)° C W^. Let X G {W^f, then D\Y){X) = 0, for ah Y eW, 
which is an immediate consequence of the definitions. But, also by definition, this 
implies that X G W^. ■ 
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Dirac structures on manifolds. We will give the definition and some basic 
properties of Dirac manifolds, following Courant [1990] and Yoshimura and Marsden 
[2006a, b] , using the notation of the latter. 

A Dirac structure D on a manifold M is a subbundle of the Whitney sum 
D C TM e T*M such that for each x G M, Dr, C T^M T*M is a Dirac structure 
on the vector space T^M. A Dirac manifold is a manifold with a Dirac structure 
on it. From Proposition 4.3 we deduce that a Dirac structure D on M yields a 
distribution -Edx C TxM whose dimension is not necessarily constant, carrying a 
presymplectic form ujd{x) ■ Edx x Edx — )■ M, for all x £ M. We can also deduce the 
following theorem, whose detailed proof appears in Yoshimura and Marsden [2006a] . 

Theorem 4.5. Let M be a manifold and let u be a 2-form on M . Given a distri- 
bution E on M, define the skew- symmetric bilinear form oje on E by restricting lo 
to E X E. For each x £ M let 

Du^EX = {{vx. ocx) e TxM e T*M I vx S E{x) 

and ax{wx) = u}e{x){vx,Wx) for all Wx G E{x)} . 

Then D^^ C TM ®T*M is a Dirac structure on AI . It is the only Dirac structure 
D on M satisfying E{x) = Edx and oje{x) = ojoix), for all x G M . 

Using Lemma 4.2 and Proposition 4.4 we can easily deduce the following propo- 
sition. 

Proposition 4.6. Let D be a Dirac structure on M and let Wx be a subspace of 
Edx for each x £ M ; then, with a self-explanatory notation, the following equalities 
hold for each x G M: 

w'x = (wf )°. 

A Dirac structure D on M is called integrable if the condition 

{Lx^a2,X3) + (Lx^asjXi) + {Lx-iai,X2) = 

is satisfied for all pairs of vector fields and 1-forms {Xi, ai), {X2, 02), (Xs, as) that 
take values in D and where Lx denotes the Lie derivative along the vector field X 
on M. This definition encompasses the notion of closedness for presymplectic forms 
and Jacobi identity for brackets. The following fundamental theorem was proven in 
Courant [1990]. 

Theorem 4.7. Let D be an integrable Dirac structure on a manifold M . Then the 
distribution E^ is involutive. If, moreover, the hypotheses of the Stefan- Sussmann 
theorem (Sussmann [1973]) are satisfied, for each x £ AI there exists a uniquely 
determined embedded submanifold S of M such that x £ S and TyS = Eoy for 
all y £ S. In other words, S is an integral submanifold of Ee,. Each integral 
submanifold S carries a presymplectic form ujd,s defined by u}£),six) = ^oix), for 
each X £ S. 



4.1 Dirac dynamical systems 



35 



4.1 Dirac dynamical systems 

Lagrangian and Hamiltonian mechanics has been developed in interaction with dif- 
ferential, symplectic, and Poisson geometry; a few references are Abraham and 
Marsden [1978]; Arnol'd [1989]; de Leon and Rodrigues [1989]; Cordero, Dodson, 
and de Leon [1989]; Marsden and Ratiu [1994]; Cendra, Marsden, and Ratiu [2001b]; 
Holm [2011a, b]. Dirac dynamical systems in the integrable case represent a synthesis 
and a generalization of both. 

Nonholonomic mechanics represents a generalization of Lagrangian and Hamilto- 
nian mechanics and is a long-standing branch of mechanics, engineering and mathe- 
matics. Some references and historical accounts on the subject are Neimark and Fu- 
faev [1972]; van der Schaft and Maschke [1994]; Zenkov, Bloch, and Marsden [1998]; 
Cendra, Marsden, and Ratiu [2001a]; Borisov and Mamaev [2002]; Cortes [2002]; 
Bloch [2003]; Cendra, Ibort, de Leon, and Martin de Diego [2004]; Cortes, de Leon, 
Marrero, and Martinez [2009]; Grabowski, de Leon, Marrero, and Martin de Diego 
[2009]; Balseiro, de Leon, Marrero, and Martin de Diego [2009]; Holm [2011a,b]. 
Some references that are more closely related to this paper are Yoshimura and 
Marsden [2006a,b]. 

Dirac dynamical systems (1.1) in the not necessarily integrable case may be 
viewed as a synthesis and a generalization of nonholonomic mechanics from the 
Lagrangian and the Hamiltonian points of view. They can be written equivalently 
as a collection of systems of the type 

i±u}D{x) = d£{x)\En. (4.1) 

A related approach has been studied in van der Schaft [1987, 1998]; van der Schaft 
and Maschke [1994]; Maschke, van der Schaft, and Breedveld [1995]; Blankenstein 
and Ratiu [2004]. As it was shown in Yoshimura and Marsden [2006a, b], one can 
write nicely the equations of nonholonomic mechanics, on the Lagrangian side, using 
the Dirac differential. We will show how this is related to system (1.1). It was 
also shown in Yoshimura and Marsden [2006a, b] and references therein how to use 
Dirac structures in LC circuit theory, on the Lagrangian side. On the Hamiltonian 
side, Poisson brackets for LC circuits were written in Maschke, van der Schaft, and 
Breedveld [1992, 1995], see also Maschke and van der Schaft [1998]. In Bloch and 
Crouch [1999] the Hamiltonian structure for nonlinear LC circuits in the framework 
of Dirac structures is investigated and simple and effective formulas are described. 

In Cendra, Marsden, Ratiu, and Yoshimura [2011] a further unification, includ- 
ing reduction, is presented, which is consistent with mechanics on Lie algebroids 
(Weinstein [1996]; Cortes, de Leon, Marrero, Martin de Diego, and Martinez [2006]; 
Martinez [2007]; Iglesias, Marrero, Martin de Diego, Martinez, and Padron [2007]; 
Martinez [2008]; Cortes, de Leon, Marrero, and Martinez [2009]). 
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5 The constraint algorithm for Dirac dynamical sys- 
tems (CAD) 

Some results in this section, in particular the CAD and examples of nonholonomic 
mechanics, are proven for not necessarily integrable Dirac structures. 

First, we shall briefly consider the case of an integrable Dirac structure. This 
includes the case of a constant Dirac structure on a vector space, that is, one that 
is invariant under translations, which includes the one used for LC circuits. Later 
on we will consider the general case. 

The case of an integrable Dirac manifold. Assume that D is an integrable 
Dirac structure on a manifold M. Each solution curve x{t) of the system (1.1) must 
satisfy {x{t),x{t)) G E£){x{t)) for all t, which implies, according to Theorem 4.7 
and equation (4.1), that it must be a solution curve to the equation 

uJD,six){x, ■ ) = i*sd£{x) (5.1) 

on a presymplectic leaf S, where is - 5" — >■ M is the inclusion, which can be solved 
using the Gotay-Nester algorithm. Such a procedure to solve (1.1) might be useful 
in those cases where the presymplectic leaves {S,u}d,s) can be found easily. This 
occurs for instance if M is a vector space and D is a. constant (i.e. translation- 
invariant) Dirac structure, as we will show next. However, this procedure has the 
drawback that in order to find a solution of (1.1) for a given initial condition, one 
must first find the leaf S containing that initial condition and then solve (5.1). For 
an initial condition on a different leaf, one has to repeat the constraint algorithm 
for a different corresponding equation (5.1). Because of this, even in these simple 
cases, working directly with the Dirac structure D, using the constraint algorithm 
to be developed in this section, rather than with the associated presymplectic form 
on a presymplectic leaf, is not only possible but also convenient, since this leads to 
obtaining a single equation on a final foliated constraint submanifold, as we will see. 

The case of a constant Dirac structure. Let y be a vector space and 
D V (B V* a given linear Dirac structure. Then we have the presymplectic form 
ujj^ on Ej^ and the associated linear map w ~ : Ej^ — )• E"- . We consider the Dirac 

structure D C TV (BT*V on the manifold V defined as Dx = {x, D), where we have 
used the natural identification TV ®T*V = Vx {V ®V*). This Dirac structure is 
integrable and constant, that is, invariant under translations in a natural sense, as 
we will show next. 

For each x €z V the presymplectic leaf Sx containing x (in the sense of Theorem 
4.7) is X + E^ C V. For each x £ V, uj^ induces the constant presymplectic form 
^^D,s £ ^'^ {x + E^) given by 

^D,s{x + y) {{x + y,X),{x + y, Y)) = Z^{X, Y), 

where y G E^, so x + y represents any point in the symplectic leaf x + Ej^, or, 
equivalently, 

io'^^six + y){x + y, X){x + y,Y)= w|(X)(F). 
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Consider the system 

{z,i)®{z,d£{z))£D. (5.2) 

This system, for a given initial condition z{0) = x, is equivalent to the following 
equation on the presymplectic leaf S that contains x 

u;D,siz)iz,-) = d£iz)\EDz. (5.3) 

Keeping x fixed and writing z = x + y, the system (5.2) on the presymplectic leaf 
S = X + becomes 

y®d£{x + y) G D. 

Therefore equation (5.3) is equivalent to the following equation on the subspace 

i:;^{y,-) = d£{x + y)\Ej^, (5.4) 

with initial condition y(0) = 0. Equation (5.4) can be solved by the algorithm 
described in Gotay, Nester, and Hinds [1978] and Gotay and Nester [1979], and 
sketched in section 2.2. Notice that, since the presymplectic form on the vector 
space Ej^ determines naturally a translation-invariant Dirac structure on the same 
space E^ considered as a manifold, equation (5.4) is also a Dirac dynamical system 
in the same way (5.2) is, but with the Dirac structure given by the presymplectic 
form ujj^ considered as a constant form on the manifold x + Ej^ . Because of this and 
also because dimE^ < dimV, equation (5.4) is, in essence, simpler than system 
(5.2). 

The case of a general Dirac manifold. Now let D he a Dirac structure on M 
that needs not be integrable. In order to explain our algorithm for Dirac manifolds 
we need the following auxiliary result, involving a given subspace Wx of E^x, which 
is easy to prove using results from section 4. 

Lemma 5.1. For each x £ M we have the following equivalent conditions, where 
Wx is a given subspace of Edx- 

(i) There exists (x,x) G Wx such that (1.1) is satisfied. 

(a) There exists {x,x) G Wx such that ujDix){x, ) = dS{x)\E£)x. 

(lii) ds{x)\EDx e wlr°. 

(iv) {d£{x)\EDx.W^^) = {Q}. 

(v) There exists {x,x) G Wx such that d£{x) £ D^{x,x). 

(vi) d£{x) G D\Wx). 

(vii) {d£{x),W^) = {0}. 
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Now we can describe the constraint algorithm for Dirac dynamical systems, 
called CAD. Following the same idea of the Gotay-Nester algorithm described in 
section 2.2 we should construct a sequence of constraint submanifolds. 

To define the first constraint submanifold Mi we may use either one of the 
equivalent conditions of Lemma 5.1, with Wx = Wqx = E^ix- We want to emphasize 
the role of the two equivalent conditions (iv) and (vii), as they represent a formal 
analogy between the CAD and the Gotay-Nester algorithm. Of course, the Gotay- 
Nester algorithm, by definition, corresponds to the case Ed = TM. 

Define 



Mi = {xeM \ {d£{x)\EDx,E'Z) = {0}} 
= {x(^M\{d£{x),E'j^\ = {Q}}. 



D^ .... (5-5) 



Let us assume that Mi is a submanifold. Then we define the second constraint 
submanifold M2 C Mi by either of the following equivalent conditions, in agreement 
with (iv) and (vii) of Lemma 5.1, with Wx = Wix = Eqx^ TxMi, 

M2 = {x G Ml I {d£{x)\EDx, {WixT^} = {0}} 
= {x G Ml I (d£{x), (Wixf) = {0}}. 



More generally we define recursively M^._|_i for k = 1, 2, . . ., by either of the conditions 

Mfc+i = {xGMk\ {dS{x)\EDx, iWkxf^) = {0}} 

= {xGMk\ (d£{x), (Wkxf) = {0}} ^^'^^ 



with Wkx = Eux^TxMk. The algorithm stops and the final constraint submanifold 
Mc is determined by the condition Mc+i = Mg. The solutions {x,x) are in Wcx = 
Edx n TxMc. 

We have proven that solution curves of (1.1) are exactly solution curves of 

{x,x)®d£[x) eDx (5.7) 
{x,x)eWcx. (5.8) 

Remark 5.2. In some examples, it is sometimes easier to write the constraint sub- 
manifolds using condition (vi) in Lemma 5.1, that is, 

Mfc+i = {x E Mfe I d£{x) e {Wkx)} 

with Wkx = Edx n TxMk, A; = 0, 1, ... , where Mq = M. 

Remark 5.3. Formula (5.6) has a special meaning in the case of an integrable Dirac 
structure. In fact, let S be an integral leaf of the distribution E^; then by applying 
the Gotay-Nester algorithm, encoded in the recursion formula (2.5), to the system 
(5.1) one obtains a sequence of secondary constraints S = Sq ^ Si ^ ■ ■ ■ ^ Sc given 
by the recursion formula 

Sk+i = {x e 5fc I {esd£{x), {TxSkT^'^) = {0}}. (5.9) 
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But it is clear that equation (5.9) coincides with equation (5.6) since T^S = Edx, 
TxSk = Wkx and the presymplectic form ujd^s on S is defined by ujd,s{x) ^ ojd{x) 
for all X € S. We are assuming regularity conditions that ensure that the Gotay- 
Nester algorithm applied for each S stops after a number of steps c which does not 
depend on S and which is at the same time the number of steps after which the CAD 
stops. As a conclusion, the final constraint submanifold Mc of the CAD is foliated 
by leaves Sc, where S varies on the set of integral leaves of the distribution Ed. We 
may say that in the case of an integrable Dirac structure D the CAD is equivalent 
to a collection of Gotay-Nester algorithms, one for each leaf S of the distribution 
Ed- The final equation given by (5.7)-(5.8) becomes 

(x,x) e d£{x) G Dx 

{x, ±) G TxSc, 

which is equivalent to the collection of equations 

uJD,s{x){x, ■ ) = i*gd£{x) 
(x, i) G TxSc- 

In section 7 we will extend the Dirac theory of constraints. For that purpose, we 
will use an embedding of M in a symplectic manifold P such that the presymplectic 
leaves S of are presymplectic submanifolds of P. The submanifold M plays the role 
of a primary foliated constraint submanifold. The case in which there is only one 
leaf gives the Dirac theory. 

Solving the equation. The constraint algorithm CAD gives a method to solve 
the IDE (1.1) which generalizes the Gotay-Nester method. Assume that the final 
constraint submanifold Mc has been determined, and consider, for each x G Mc, the 
affine space 

S^^) ■= {{x,±) G TxMc I (4.1) is satisfied}, 

which is nonempty if Mc is nonempty, a condition that will be assumed from now 
on. Let d^'^\x) be the dimension of si^\ 

The following theorem is one of the ingredients of our main results, and gener- 
alizes the Gotay-Nester algorithm for the case of Dirac dynamical systems (5.10) 
rather than Gotay-Nester systems (1.2). Its proof is not difficult, using the previous 
lemma, and is left to the reader. 

Theorem 5.4. Let M be a given manifold, D a given Dirac structure on M and E 
a given energy function on M , and consider the Dirac dynamical system 

{x,x)®d£{x) e Dx. (5.10) 

Assume that for each k = 1, ■ ■ ■ , the subset M^ of M defined recursively by the 
formulas (5.5)-(5.6) is a submanifold, called the /c-constraint submanifold. The 
decreasing sequence M^ stops, say Mc = Mc+i (which implies Ale = Mc+p, for all 
p and call Mc the final constraint submanifold. Then the following hold: 
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(a) For each x S Mc, there exists {x,x) € Wcx = H T^Mc such that (5.10) 
is satisfied. The Dirac dynamical system (5.10) is equivalent to the equation 

ujd{x){x , ) = d£{x)\EDx, {x,x) G Wcx, (5-11) 

that is, both equations have the same solution curves x{t) G Mc- 

(b) For each x £ Mc, d^'^\x) equals the dimension ofkeiujDix) CiTxMc. 

(c) If d^'^\x) is a locally constant function of x on Mc then S^^^ = U^eAfc '^^'^^ 
is an affine bundle with base Mc. Each section X of 5^^^ is a vector field on Mc 
having the property that X(x)(Bd£{x) G Dx, for all x G Mc. Solution curves to such 
vector fields are solutions to the Dirac dynamical system (1.1). More generally, one 
can choose arbitrarily a time- dependent section Xt, then solution curves of Xt will 
be also solutions to (1.1) and those are the only solutions of (1.1). Solution curves 
to (1.1) are unique for any given initial condition if and only if d^^\x) = 0, for all 
X G Mc. 

Local representation for 5^^^^. One can find a local representation for S^^^ by 
just choosing a local parametrization of Mc. Let x = x{zi, . . . , Zr) G U, where U is 
an open set and r is the dimension of Mc, be such a local parametrization. Then 
substitute this expression for x in (1.1) to obtain an IDE m z = (zi, . . . , Zr), namely 

{x{z),D,x{z) ■i)®dS {x{z)) G D^(,y (5.12) 

The local representation of S^'^'' is given by the trivial affine bundle 

{{z,z) I {z,z) satisfies (5.12)}. 

Some results concerning uniqueness of solution. Under Assumption 3.18 
(b), we can prove the following lemma. 

Lemma 5.5. (A) Existence and uniqueness of a solution curve x{t) of (5.11) for 
any initial condition x{Q) G Mc, and therefore also of the Dirac dynamical system 
(5.10), is equivalent to any of the conditions 

(i) k.Qiujj:){x) n Wcx = kera;£)(x) n TxMc = {0}, for each x G Mc, 

(ii) Wcx n = {0}' for each x G Mc, 
(Hi) Wcx n E%^ = {0}, for each x G Mc. 

(B) If ooo\Wcx, or, equivalently, Q.\Wcx, is symplectic for each x G Mc then 
there is existence and uniqueness of solution x{t) of (5.11) for any initial condition 
x{0) G Mc, and therefore also of the Dirac dynamical system (5.10). One also has 
that uJo\Wcx is symplectic for each x G Mc if and only if any one of the following 
conditions is satisfied: 

(i)WcxnWiff = {0}, 
(zt)Wcxnw^, = {0}. 
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Proof. We first recall the argument from Theorem 5.4. Uniqueness of a solution 
x{t) of equation (5.11) for any given initial condition xq G Mc (which, as we know 
from the CAD, must satisfy x{t) £ Mc and x(t) £ W^x^t) for all t) holds if and only 
if for each x G Mc there is a uniquely determined Vx G Wcx such that 

UJDix){Vx, ) = d£\EDx, 

in other words, if equation (5.11) defines a vector field on Mc- In fact, if this is the 
case, by the general theory of ODE on manifolds we have existence and uniqueness of 
solution. We also recall that under Assumption 3.18 (b), the equation (5.11) defines 
a family of vector fields on Mc, which defines an affine distribution of constant rank 
whose space at the point x G Mc is the affine space of all solutions Vx as indicated 
above. Using this we can deduce that uniqueness of solution x{t) of (5.11) for any 
initial condition x(0) G Mc, is equivalent to the affine distribution defined above 
having dimension 0. 

Now we shall prove the equivalence between uniqueness of solution and (A). Let 
Vx G Wcx satisfying equation (5.11), that is, 

0JDix)ivx, ) = d£\EDx- 

Let Wx 7^ be such that Wx G kerw£)(x) n Wcx- Then Vx + Wx also satisfies (5.11), 
which shows that the affine distribution described above has dimension greater than 
0. Using this, the proof of the equivalence between uniqueness of solution and (A)(i) 
follows easily. The rest of the proof of the equivalence with (A) follows from the 
fact that 'k.ei u}d{x) = = E^^ n E^x for each x G Mc, which can be proved 
directly using the definitions. 

Note that the equivalence between (i), (ii), (in) oi (A) holds for any subspace 

Wcx ^ Edx- 

Now we shall prove (B). First of all, if Vx G Wcx satisfies (5.11) then, since 
Wcx C E£)x, it clearly also satisfies 

{ujd{x)\Wcx){vx, ) = d£\Wcx. (5.13) 

We can conclude that if u:d{x)\Wcx-, which as we know coincides with Q.{x)\Wcx, is 
symplectic then equation (5.11) defines a vector field Vx on Mc, which is, in fact, 
given by equation (5.13). Then existence and uniqueness of solution x{t) of (5.11) 
for any initial condition x(0) G Mc is guaranteed. Finally, the equivalence between 
symplecticity oioj£){x)\Wcx = f^(2;)|VFcx and (B)(i) and (B)(ii) is easy to prove using 
basic linear symplectic geometry. In fact, is is easy to prove that symplecticity of 
ijjDix)\Wcx is equivalent to = {0}. Using Lemma 3.4 we can deduce that 

W^x" = Wcx n Wc^ from which we obtain Wcx n Wg = Wcx n W^x"" ■ ■ 

With the method just described, one can deal with many examples of interest, 
such as nonholonomic systems and circuits, provided that one chooses the manifold 
M and the Dirac structure D properly. We show in the next section how a nonholo- 
nomic system given by a distribution A C TQ on the configuration space Q can be 
described by a Dirac dynamical system on the Pontryagin bundle M = TQ ®T*Q 
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and how one can apply the constraint algorithm CAD to this example, although 
we will not perform a detailed calculation of the sequence of constraint submani- 
folds Mfc. We also show how LC circuits can be treated with the same formalism 
as nonholonomic systems. The main point for doing this is, again, to choose the 
manifold M as being the Pontryagin bundle TQ (B T*Q, where, this time, Q is 
the charge space, and a canonically constructed Dirac structure -Da on M, where, 
this time, A represents Kirchhoff 's Current Law. We also show how this approach 
using TQqT*Q is related to the approach used in Yoshimura and Marsden [2006a, b]. 



6 Examples 

In this section we deal with two examples, namely, nonholonomic systems, and LC 
circuits, showing that Dirac dynamical systems give a unified treatment for them. 
We will perform the detailed calculation of the constraint submanifolds Mj. for the 
case of LC circuits only. 

Nonholonomic systems. Recall that a nonholonomic system is given by a con- 
figuration space Q, a distribution A C TQ, called the nonholonomic constraint, and 
a Lagrangian C: TQ — )• M. Equations of motion are given by Lagrange-d'Alembert's 
principle. 

Inspired by the Hamilton-Poincare principle given in Cendra, Marsden, Pekarsky, 
and Ratiu [2003a], we can write a convenient equivalent form of the Lagrange- 
d'Alembert principle as 

6 / {pq - £{q,v,p))dt = 0, 

Jto 

where £: TQ ® T*Q — )• M is defined by £{q,v,p) = pv — C{q,v), and with the 
restriction on variations 6q G A, 5q{ti) = for z = 0, 1, along with the kinematic 
restriction v £ A. The resulting equations are 



AO 

oq 


(6.1) 


q = V 


(6.2) 


dv 


(6.3) 


V e A. 


(6.4) 



We are going to show that equations (6.1)^(6.4) can be written in the form (1.1). 
For this purpose we must construct an appropriate Dirac structure associated to 
the nonholonomic constraint. Inspired by several results in Yoshimura and Marsden 
[2006a] and by the Hamilton-Poincare's point of view we define a Dirac structure 
Da ^ TM eT*M on M = TQ eT*Q associated to a given distribution A CTQ 
on a manifold Q by the local expression 

DA{q,v,p) = {{q,v,p,q,v,p,a,'y,P) \ q G A{q), a + pe A°{q), /3 = g, 7 = 0}. 
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Note. We shall accept both equivalent notations 

{q,v,p,q,v,p,a,j,^) = {q,v,p,q,v,p) {q,v,p,a,j, P), 

for an element of TM T*M. 

By checking that diinD/\{q,v,p) = SdimQ, that is dim D/\{q,v,p) = dimM, 
and that 

{{{q,v,p,qi,vi,pi) e {q,v,p,ai,ji,Pi), {q,v,p, q2,V2,P2) © (g, a2 , 72, /32)» = 
for all 

{Q,v,p,qi,vi,Pi) © (g, v,p,ai,7i,/3i) G D/^{q,v,p), 
{Q,v,p,q2,V2,P2) © {q,v,p,a2,-t2,l32) G D/^{q,v,p), 

we can conclude using Lemma 4.1 that -Da is a Dirac structure on M. We should 
now prove that Z)^ is well defined globally, in other words, that it does not depend 
on the choice of a local chart. Let f : TQ r*Q ^ Q and vf : TQ T*Q ^ T*Q be 
the natural maps that in local coordinates are given by f{q, v,p) = q and 7f(g, v,p) = 
{q,p). For a given distribution A C TQ consider the distribution A = (Tr)~^(A) 
and also the 2-form O = Tt*Ljj, on the manifold TQ ®T*Q, where oj is the canonical 
2-form on T*Q. We have the local expressions A = {{q,v,p,q,v,p) \ q G A} and 
Q{q,v,p) = dq A dp. Now we can apply Theorem 4.5 replacing M by TQ © T*Q, 
£^ by A and lv hy O and then we can easily check that the Dirac structure Di^^ 
coincides with D/^. In other words, by using Theorem 4.5 we have obtained a 
coordinate- independent description of D/^ given in terms of f, vf and O, which 
proves in particular that it is well defined globally. 

It is straightforward to check that the condition 

{x,±)ed£{x) e Da, (6.5) 
where x = {q,v,p), is equivalent to 

q = V (6-7) 

4 6 A, (6-9) 

which is clearly equivalent to equations (6.1)"(6.4). 

Now that we have written the equations of motion as a Dirac dynamical system, 
we can proceed to apply the CAD. First, we can easily prove the following formulas 
using the definitions and Proposition 4.4: 

= {(g,z;,p,a,7,/3) |7 = 0,/3G A}, (6.11) 
={(g,z;,j),g,i>,p) \q = 0,peA°}. (6.12) 
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Then we have 

Ml = [{q,v,p) I (^d£{q,v,p),E^2('l^v,p)) = {0}} 



\^{q,v,p) d£{q,v,p) e E^^^{q,v,p)^ (513) 
{q,v,p) 



— = 0, VGA 
av 

We could continue applying the algorithm as explained in general in section 5, 
and we would obtain specific formulas for M2,M3, etc. 

Relationship with implicit Lagrangian systems. In Yoshimura and Mars- 
den [2006a] the description of a nonholonomic system as an implicit Lagrangian sys- 
tem was introduced and equations of motion were shown to be a partial vector field 
written in terms of the Dirac differential. There is a close and simple relationship 
between this approach and the one of the present paper, which we shall explain 
next. 

First, we should recall the notion of an implicit Lagrangian system. Let C : TQ — )• 
M be a given Lagrangian; then dC is a 1-form on TQ, dC : TQ — )• T*TQ, which is 
locally expressed by 

/ dC dC\ 
\ oq ov J 

Define a differential operator 2), called the Dirac differential of C, by 

D£ = 7Q o dC, 

where 7q is the diffeomorphism defined by 

7Q = 0^ o : T*TQ ^ T*T*Q. 

Here Vt is the canonical 2-form on T*Q and kq: TT*Q — )• T*TQ is the canonical 
isomorphism which is given in a local chart by KQ{q,p, 6q, dp) = {q,6q,6p,p). We 
have the local expression 

/ dC dC 
'i^C= [q, — , -^,v 
\ ov oq 

We must now recall the definition and properties of the Dirac structure on 
the manifold T*Q, where A is a given distribution on Q, studied in propositions 
5.1 and 5.2 of Yoshimura and Marsden [2006a]. Let tt: T*Q — )• Q be the canonical 
projection. Then At*q is defined by 

AT.Q = (T7r)-i(A). 

The Dirac structure -Da is the one given by Theorem 4.5, with M = T*Q, lo = 0, 
the canonical symplectic form on T*Q and E = Aj-g, and it is described by 

DA{q,p) = {v®ae ^,,p)T*Q Tl^^^.^T*Q \ v G At.q((Z,p) 

and a{w) = il{q,p){v,w) for all w G AT*Q{q,p)}- (6-14) 

Now we will define the notion of an implicit Lagrangian system. 
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Definition 6.1. Let C: TQ — )■ R &e a given Lagrangian (possibly degenerate) and 
let A C TQ be a given regular constraint distribution on a configuration manifold 
Q. Denote by Da the induced Dirac structure on T*Q that is given by the equation 
(6.14) and write T):TQ^ T*T*Q for the Dirac differential of C. Let FL(A) C 
T*Q be the image of A under the Legendre transformation. An implicit Lagrangian 
system is a triple (C, A, X) where X represents a vector field defined at points of 
FL(A) C T*Q, together with the condition 

X®T)C£ Da- 

In other words, for each point (q, u ) G A let {q,p) = ^L{q, v) and then, by definition, 
X{q,v,p) £ T(^q,p) iT*Q) must satisfy 

X{q,v,p)®T)C{q,v) e DA{q,p). (6.15) 

A solution curve of an implicit Lagrangian system (C, A, X) is a curve {q{t), v{t)) £ 
A, ti < t < t2, which is an integral curve of X where {q{t),p{t)) = ¥L{q{t),v{t)). 
The following proposition is essentially Proposition 6.3 of Yoshimura and Marsden 
[2006a]. 

Proposition 6.2. The condition X © DC £ defining an implicit Lagrangian 
system is given locally by the equalities 

P=f (0.10) 

qeA (6.17) 
v = q (6.18) 

p-|^GA°. (6.19) 

It is clear that equations (6.16)-(6.19) are equivalent to (6.1)~(6.4) and also 
to (6.6)-(6.9). This leads immediately to a precise link between the approach to 
nonholonomic systems given in Yoshimura and Marsden [2006a] and the one in the 
present paper, which is given in the next proposition whose proof in local coordinates 
is easy and will be omitted. Recall that Mi is the first constraint manifold given by 
(6.13). 

Proposition 6.3. Let C: TQ — t- M &e a given Lagrangian, A C TQ a given distri- 
bution and let E: TQ © T*Q — t- M 6e given by £{q,v,p) = pv — C{q,v). Then the 
following assertions hold. 

(i) FL(A) = 7f(Mi). 

(ii) Let (q, v,p, q, v,p) = X{q, v,p) be a solution to (6.5); then in particular one has 
(g, v,p) G Ml, and let X{q, v,p) = {q,p, q,p), that is, X{q, v,p) = TTtX{q, v,p). 
Then X{q,v,p) is a solution to (6.15). 

(iii) A curve {q{t),v{t),p(t)) is a solution curve of (6.5) if and only if {q{t),v{t)) 
is a solution curve of the implicit Lagrangian system given by (6.15). 
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LC circuits. This is a case of a constant Dirac structure of tlie type explained 
in section 5. Our approacii is closely related and equivalent to the one described 
in Yoshimura and Marsden [2006a]. There are many relevant references where the 
structure of Kirchhoff 's laws has been studied from different points of view; some 
of them emphasize the geometry behind the equations, see Smale [1972]; Chua and 
McPherson [1974]; Maschke, van der Schaft, and Breedveld [1995]; Moreau and 
Aeyels [2004]. 

We first briefly recall the description of LC circuits given in Yoshimura and 
Marsden [2006a]. Let E he a vector space representing the charge space; then TE 
is the current space and V = T*E is the flux linkage space. There is a constant dis- 
tribution A C TE (that is, A is invariant under translations), which represents the 
Kirchhoff 's Current Law (KCL). The subbundle A° <Z T*E represents the Kirch- 
hoff 's Voltage Law (KVL). One has a Dirac structure on the cotangent bundle T*E, 
Da^TV ®T*V given by 

Da = {{q,P,q,p) © iq,p,ag,ap) £ TV ®T*V \ q £ A, Op = q, ag + p £ A°}, 

which clearly does not depend on {q,p), in other words it is a constant Dirac struc- 
ture, so for each base point {q,p) we have 

DA{q,p) = {q,p,D), 

where 

D = {{q,p,a,l3)eV®V* \qeA, p = q,a+p£A''}. 

The dynamics of the system is given by a Lagrangian C : TE — t- M. This Lagrangian 
is given by a quadratic form on E x E representing the difference of the energies in 
the inductors and the capacitors, say, 

1 " 1 " 1 

1=1 i=l * 

where n is the number of branches of the circuit. Of course, some of the terms in 
the previous sum may be zero, corresponding to the absence of an inductor or a 
capacitor in the corresponding branch. 

One approach would be to use the time evolution of the circuit in terms of the 
Dirac differential, given by equation (6.15). This shows that Dirac structures provide 
a unified treatment for nonholonomic systems and LC circuits. In other words the 
equation is 

{¥C{q,v),q,p) © {q,p,^C{q,v)) £ DAiq,p), 

with {q,p) = ¥C{q,v) and {q,v) £ A, where DC represents the Dirac differential of 
£. 

However, we will follow the philosophy of this paper and choose to work on 
TE (BT*E rather than working with the Dirac differential. 
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Applying the CAD to LC circuits. It should be now clear how we can 
apply the constraint algorithm CAD developed in the present paper to deal with 
LC circuits exactly as we did with the case of nonholonomic systems. More pre- 
cisely, it should now be clear that an LC circuit can be described by the methods 
described in the first part of the paragraph Nonholonomic Systems (beginning of 
section 6), by taking Q = E, M = TE®T*E and defining Z)a and 8 as indicated in 
that paragraph. As we have already said, in the case of circuits the Dirac structure 
-Da is constant and therefore integrable. However, we prefer not to work with the 
system restricted to a presymplectic leaf, as explained before, but to apply directly 
the algorithm for a general Dirac structure, as we will explain next. 

Define the linear maps (p: E ^ E* and ip: E ^ E* hy 

dC 

(p{v) = — = {LiVi, . . . , LnVn) , 

dC 

i^iq) = ^ = ~ {Qi/Ci, . . ■,qn/Cn) ■ 

If there is no capacitor in a circuit branch, the corresponding component of ip{q) will 
be zero (infinite capacitance), while zero capacitance is ruled out since it would rep- 
resent an electrically open circuit branch. The physically relevant case corresponds 
to < Lj < oo, < Cj < oo, for i = 1, . . . ,n. However, the constraint algorithm 
applies to the general case where capacitances and inductances can be negative. 

The evolution equations (6.6)~(6.9) for a general nonholonomic system become 

p-ip{q)£A° (6.20) 

q = v (6.21) 

p = ip{v) (6.22) 

qe A. (6.23) 

Now we will apply the CAD. First, note that 

d£{q, v,p) = ip{q)dq + {p — ip{v))dv + vdp. 

The first constraint submanifold is calculated as follows, taking into account the 
general expression (6.11), 

Ml = {[q,v,p) £ M I d£{q,v,p) £ E^jyjq,v,p)} 
= {{q,v,p) I P = ^{v), V G A}. 
Now we simply continue applying the algorithm. Let 

Wi=TMinEQ^ ={{q,v,p,q,v,p) \ {q,v,p) e Mi, q € A, p = ^{v), i) € A}. 
Then l)^(VFi), denoted simply as W^, is 

Wt = {{q,v,p,a,j,l3) \ {q,v,p) £ Mi, a £ (^(A) + A°, 7 = 0, {3 £ A}, 
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and therefore 

M2 = {{q^v^p) G Ml I d8{q,v,p) S Wl{q,v,p)] 

= {{q.v.p) I q G V'^Hv'CA) + A°), p = ip{v), v E A}. (6.24) 

In the same way we can calculate M3. In fact, 

W2=TM2nEjj^={{q,v,p,q,v,p) \ {q,v,p) e M2, q e A,q e iJ~\ip{A) + A°), 

P = 'Piv), i) G A}, 

= {(g,T;,p,a,7,/3) | {q,v,p) e M2, a e V9(A)+A°, 7 = 0, /3 E Ar\^-\^{A)+A°)] , 
then 

= {{q,v,p) G M2 I d8{q,v,p) S 1^2^(5, 

= {(9, v,p)\qe r\v{A) + A°), p = ^{v), v€An ^-HifiA) + A°)} 
= {{q,v,p) \ q e ij'\ipiA) + A°), p = ifiv), V £ Ai}, 

where we have called 

Ai = A n ij-\ip{A) + A°) c A. 

Now we shall calculate M4. We will show later that M4 = M3 for all LC circuits with 
either a positive inductance or a non-infinite positive capacitance on every branch 
(Lemma 6.5 and Theorem 6.6). First we compute 

W3 = TM3 n E^^ 

= {(.q,v,p,q,v,p) I {q,v,p) G M3, qe A,qe il)~^{ip{A) + A°), p = ip{v), v G Ai} 
= {{q,v,p,q,v,p) I {q,v,p) G M3, q e Ai, p = ip{v), v G Ai}. 

Thus, one finds that 

= {{q,v,p,a,-f,P) I {q,v,p) e M3, a £ (^(Ai) + A°, 7 = 0, /? G Ai}, 
and hence 

M4 = {iq,v,p) G M3 I d£{q,v,p) G W3^(g,w,p)} 

= {(g,f,p) G M3 I g G V'-i(^(Ai) + A°), p = ip{v), v G Ai} 
= {iq,v,p)\q£ V~i(vp(Ai) +A°),p = ^{v), v G Ai}. 

Similarly we calculate M5 as follows. 

W4 = TM4 n -E^,^ 

= {{q,v,p,q,v,p) I (g,u,p) G M4, g G A,g G V~^(v?(Ai) + A°), p = ip{v), v G Ai} 
= {{q,v,pA,v,p) I (9,^^,P) e M4, g G AnV'-H(^(Ai) + A°), p = 99(7)), i> G Ai} 



6 Examples 49 



Let A2 = A n t/' ^{(f{Ai) + A°); then A2 C Ai and we can rewrite W4 as 
W4, = {{q,v,p,q,v,p) I {q,v,p) G M4, q e A2, p = ^{v), v G Ai}. 
Then we have 

Wl = {{q,v,p,a,-i,l3) I (Q,i;,p)GM4, a G ^(Ai) + A°, 7 = 0, ^ G A2} 

= G M4 I G Wl{q,v,p)] 

= {{q,v,p) G M4 I g G '0"^('^(Ai) + A°), p = if{v), v G A2} 

= {{q,v,p) \ q e -^-^^iAi) + A°), p = ^{v), V e A2}. 

From this, one sees how to recursively define Mk- For all /c > 1 define 

Afc = Anv^-^MAfc^O + A"), 

where Aq = A by definition. We have the following expressions for the constraint 
submanifolds Mk- 

Ml = {iq,v,p) I p = ^p{v),v G A} 
M2fc_i = {iq,v,p) I q £ ij-\ipiAk-2) + A°), p = ^{v), v G Afe_i}, k>2 
M2k = {iq,v,p) I q £ ij-\ipiAk-i) + A°), p = ^{v), v G Afe_i}, A; > 1. 

In order to solve the system we may apply the method to solve the general equa- 
tion (5.12) explained before. It is clear that the parametrization x = x(zi, . . . , Zr), 
where r is the dimension of the final constraint Mc, can be chosen to be a linear 
map and the IDE (5.12) will be a linear system. 

This system is of course still an IDE. One of the purposes of this paper is to 
write Hamilton's equations of motion, using an extension of the Dirac procedure 
extended to primary foliated constraint submanifolds, developed in section 7. For 
this purpose, we will need to consider M = TE ®T* E as being the primary foliated 
constraint embedded naturally in the symplectic manifold T*TE. 

Physical interpretation of the constraint equations for LC circuits. 

The equations defining the constraint submanifolds have an interesting interpreta- 
tion in circuit theory terms. From now on we will assume the physically meaningful 
situation where < Lj < 00, < Cj < 00, i = 1, . . . , n. We will show that in this 
case, the algorithm stops either at Mi or at M3. 

A circuit has an underlying directed graph, since each branch has a direction 
in which the current flow will be regarded as positive. A loop is a closed sequence 
of different adjacent branches. It also has a direction, which does not have to be 
compatible with the directions of the branches involved. The set of loop currents is 
a set of generators for the KCL subspace ACE'. 

We shall call a circuit branch inductive (resp. capacitive) if there is an inductor 
(resp. capacitor) present on that branch. A loop will be called inductive (resp. 
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capacitive) if it has an inductor (resp. capacitor) on at least one of its branches. A 
branch or loop that is capacitive but not inductive will be called purely capacitive. 
A branch or loop that is neither inductive nor capacitive will be called empty, since 
no capacitors or inductors are present. Thus, a non-inductive branch or loop can be 
either empty or purely capacitive, and a purely capacitive loop must have at least 
one capacitor, no inductors, and possibly some empty branches. 

The first one of the equations defining Mi is p = (p{v), which is pi = LiVi, 
i = 1, . . . ,n. The quantity LiVi is the voltage corresponding to an inductor through 
which the current is Uj. Then pi can be interpreted as a time integral of the voltage 
on branch i due to the inductor in that branch (the flux linkage of the inductor). 
The second equation, G A, is just KCL for the currents v. 

The next constraint submanifold M2 incorporates the equation ijj(q) G 97(A) + 
A°, which represents the KVL equations for the purely capacitive loops, as explained 
in the next theorem. 

Theorem 6.4. The subspace (p{A) + A° C £;* represents Kirchhoff's Voltage Law 

for non-inductive loops (KVLNI for short). That is, its elements are precisely those 
branch voltage assignments on the circuit that satisfy the subset of KVL equations 
corresponding to the non-inductive loops. The equation ipiq) G (^(A) + A° represents 
the condition that on every purely capacitive loop, the branch voltages qi / Ci of the 
corresponding capacitors satisfy KVL. We call these the KVLPC equations (KVL 
for purely capacitive loops), and they form a subset of the KVLNI equations. 

Proof. Let {ei, . . . , e„} be the basis of E where is associated to branch i of the 
circuit, and let {e^, . . . , e"} be its dual basis. Recall that the branch voltage assign- 
ments on the circuit are elements of E*, and A° represents Kirchhoff's Voltage Law 
in the sense that each element of A° is a branch voltage assignment satisfying KVL. 
The KVL equations are linear equations on the branch voltages and are therefore 
represented by elements of E. In fact, they are precisely the elements of (A°)° = A. 

Consider a maximal set of independent, non-inductive oriented loops on the 
circuit, labeled l,...,m. Each loop gives rise to a linear equation on the corre- 
sponding branch voltages according to KVL. These KVL equations are represented 
hy r]i, ... ,r]m & E, and are a subset of the full KVL equations for the circuit. That 
is, these linear equations hold on any branch voltage assignment that is compatible 
with KVL, and therefore rji,.. . ,r)m € A. Note that each rjj = rjjBi does not involve 
any inductive branches; that is, if 7^ then rjj = r]j{e^) = 0. Also, the KVLNI 
equations r]i, . . . .rjm can be seen as the KVL equations for the circuit that is ob- 
tained by removing the inductive branches from the given circuit. Let us call this 
the non-inductive subcircuit, whose charge space £'ni is a subspace of in a natural 
way. Also, E^^ C E* is generated by {e*} where i ranges over all non-inductive 
branches. 

Let us denote Ani C the KCL distribution for the non-inductive subcir- 
cuit. Note that Ani C A, since any KCL-compatible current assignment on the 
non-inductive subcircuit corresponds to a KCL-compatible current assignment on 
the original circuit for which no current flows through the inductive branches. Re- 
ciprocally, elements of A that are zero on the inductive branches can be regarded 
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as elements of Ani- Note that rji, . . . ,rjm are loop currents for the non-inductive 
subcircuit, and they form a basis of Ani- 

For any v = v^Ei G E, ip{v) = ^ LiV^e^. Then for each j = 1, . . . ,m, r]j{(p{v)) = 
^Lif*?7j(e*) = 0. This is true in particular for v e A., so rjj G {ip{A))°. Therefore 
rfj G {ip{A))° n A = {ip{A) + A°)°. In other words, each 771, ... , rjm, seen as a linear 
equation on E*, holds on ip{A) + A° C E*. 

Let us see that (p{A) + A° is precisely the vector subspace of E* defined by 
the equations r/i, . . . ,rim- Let p = p^ti G ('/'(A) + A°)° = {ip{A))° n A be another 
equation that holds on <^(A) + A°. Then = p{(p{p)) = J2^iiP^)^- Since all 
inductances are nonnegative, this implies that the components of p corresponding 
to the inductive branches must be zero. Also, /? G A, therefore p G Ani- This means 
that p is a linear combination of 771, . . . , r/m- 

The equation ijj{q) G ip{A) + A° means that r]j{ilj{q)) = for j = l,...,m. 
Only those rjj corresponding to purely capacitive loops give a condition on q. The 
remaining r]k correspond to empty loops, so r]k{i/j{q)) = amounts to the equation 
= 0. ■ 

The third constraint submanifold M3 is obtained from M2 by incorporating the 
equation i/j{v) G '/?(A) + A°, which means that the currents on the branches of the 
purely capacitive loops satisfy the same KVLPC equations as the charges. It is 
clear that it must hold when we consider the dynamics, which includes the equation 
q = V. 

Suppose that there exists at least one purely capacitive loop so that the algo- 
rithm does not stop at Mi. Then it is possible to show that M3 C M2. In order to 
do that, consider a purely capacitive loop, which by definition must involve at least 
one capacitor, and recall that all the capacitances are positive. The corresponding 
KVLPC equation is represented by rj e E, whose components are ±1 or 0. As in 
the proof of the previous theorem, rj G (A°)° = A, so it can be interpreted as a 
nonzero loop current, that is, a current that is ±1 on all the branches of that loop 
(depending on the relative orientation of the branches with respect to the loop) and 
zero on the remaining branches. Note that V'(^) G E* does not satisfy KVLPC, 
since in particular r){ip{r])) = ^{rfY /Ci = ^ l/Cj > 0, where the last sum is over 
the branches in the chosen loop. Therefore A ^ ^~^((^(A) -|- A°) and M3 C M2. 

Now we will prove the interesting fact that since inductances are greater than 
or equal to zero and capacitances are positive numbers or -|-oo the algorithm always 
stops at M3 (or earlier). We shall start by proving the following lemma. 

Lemma 6.5. (a) 61 G {(p{A) + A°) is equivalent to the condition that there 
exists (5o G A such that tp{Si) — ip{So) G A°. 

(b) The condition Ai = A2 is equivalent to the following condition: 



Condition (*) 
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For any 5 € A there exist 5q,5i S A such that the following hold: 



ip{5)-ip{5i)e^° 
V'(5i)-(/^(<^o)e A°. 



Equivalently, we can say that for each 5 € A there exists 6i S Ai such that 



(c) Condition (*) implies that if = M5 then M3 = M4. 

Proof. The proof of (a) is immediate taking into account the definitions. To prove 
(b) take any given 5i G Ai, then it satisfies the condition V'('^i) ~ ^{^o) ^ 1 
where Sq E A. By condition (*) there exists S'^ S Ai such that if {5q) — ip {S'l) G A°, 
and we can conclude that V'('^i) ~ ¥'('^'1) ^ ■> which shows that 5i G A2. The 
converse can be easily verified. 

Now we shall prove (c). Since the inclusion M4 C M3 is immediate by con- 
struction we will only prove the converse. Let {q,v,p) S M3 then, in particular, 
there exists 5 G A such that "0(9) ~ V'i^) ^ By condition (*) there exists 

61 £ Ai such that if {61) — if{5) G A°, which implies that ^p{q) — 95(^1) G A°, that 
is, q G tp'^ {ip{Ai) + A°)), from which we can deduce that M3 C M4. ■ 

The following theorem gives a sufficient condition under which condition (*) in 
the previous lemma holds, which considers the standard physically meaningful case. 
Also, we assume for simplicity that there are no empty branches. We should mention 
that it can be proven that even with empty branches, the algorithm stops at Mi or 
Ms. 

Theorem 6.6. Condition (*) in Lemma 6.5 holds if in every branch there is either 
a positive inductance or a non-infinite positive capacitance. 

Proof. Condition (*) can be reformulated as follows: 

For any 6 £ A there exists 5o,5i G A such that the following equalities are 
satisfied 



for all S[,6q G A. This formulation has the advantage that := ip + ^p is nonsingular. 

The basis Cj, i = 1, . . . ,n of the charge space E defines naturally an Euclidean 
metric by the condition (cj, ej) = 5ij and an identification E = E* by the condition 
Cj = e', i = 1, . . . , n, in particular obtains A° = A-*-. The linear map v is self-adjoint 
and positive definite while if is self-adjoint and positive semi-definite. System (6.25)- 
(6.26) can be written in the form 



ip{5)-ip{6i)eA°. 



((V^ + (^) {Si),5',)-{^{6o),5'o) = {f{6),6'o) 



(6.25) 
(6.26) 



PAOip\A{6i) = PAOip\A{6) 
Pa o i^\A{5i) - Pa o V'|A(<5o) = Pa « ^|A(<5o). 



where Pa : E ^ A is the orthogonal projection on A. 
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One can choose an orthonormal basis of A and then the system of equations 
(6.25)-(6.26) can be written in matrix form. Moreover we have an orthogonal de- 
composition A = A^^^ + A^^^ where Ai and A2 are the image and the kernel of 
Pa o V'IA. Each vector 5 G A is decomposed as = 6^^^ + 6^'^\ We can choose 
the basis in such a way that the block decomposition of the matrix representing 
Pa o '/'IA has the form 



(1,1) ^(1,2) 
(2,1) ,^(2,2) 



where ip 



(1,2) 



(2,1) 



(2,2) 



,(1,1) 



and ip^ ' ' is a diagonal matrix with positive 



eigenvalues. The map Pa o z^|A also has a block decomposition, and since it is self- 
adjoint and positive definite, the blocks u^^'^') and i^^^'^^ are symmetric and positive 
definite matrices. We obtain the following system of equations, which is equivalent 
to (6.25)-(6.26), 



(1,1) (^(1,2) 
(2,1) y,(2,2) 



6 



(1) 

.(2) 



(1,1) c^(l,2) 
(^(2,1) (^(2,2) 



1,(1,1) ^(1,2) 
1.(2,1) ^(2,2) 





\ ^1^ 1 






_ <5f ) _ 





(1,1) (^(1,2) 
(2,1) (^(2,2) 



For given 5, equation (6.27) fixes jj^^ and imposes no condition on 6\^' . Using the 
fact that if^^'^'^ and z^^^'^^ are invertible we can find and (5f ^ to satisfy (6.28). ■ 

As a conclusion, if (and only if) there are no purely capacitive loops in the 
circuit, then M2 = Mi and the final constraint submanifold is Mi, defined by the 
conditions that pi = LiVi (where some Li might be zero) and v satisfies KCL. 
Otherwise, the final constraint submanifold is M3, defined by the conditions that 
Pi = LiVi, V satisfies KCL, and qi/Ci and Vi/Ci satisfy KVLPC. Recall that qi/Ci is 
the voltage of the capacitor in branch i, and the absence of a capacitor on branch i 
means Cj = 00. The rest of the KVL equations, which involve the branch voltages 
Pi = LiVi on the inductors, will appear when considering the dynamics. In this 
sense, the KVLPC equations can be regarded as "static" KVL equations. 

Remark 6.7. The hypothesis that no capacitances or inductances are negative, be- 
sides being physically meaningful, is crucial to ensure that the algorithm stops at 
M3 (or Ml). For example, for a circuit with one inductor L, one capacitor C and a 
negative capacitor — C, all of them in parallel, it stops at M5. 



5. 



(1) 


^(2) 



(^(2,1) (^(2,2) 



(2) 



(6.27) 

5(1) 

5(2) 

(6.28) 



Geometry of the final constraint submanifold. We know that the solution 
curves will have tangent vectors in Wc = TMcCiEi)^ , where c = 1 or c = 3 depending 
on the case. 

Theorem 6.8. Assume the physically meaningful situation where < Li < 00, 
< Cj < 00, i = 1,. . . ,n. The distribution Wc C TMc is constant and therefore 
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integrable, and its integral leaves are preserved by the flow. Denote the leaf through 
by Wc, so these leaves can be written as x + Wc, where x = {q,v,p). Each integral 
leaf X + Wc is a symplectic manifold with the pullback of the presymplectic form 
LO = dq A dp, if and only if there are no empty loops. 

Proof. As mentioned before, if there are no purely capacitive loops, then c = 1, 
otherwise c = 3. For the first case, using (6.10), we have 

Wi = {{q,v,p,q,ij,p) e TEeT*E\{q,v,p,q,v,p) eTMi,qeA} 

= {{q, v,p, q, v,p) G TE © T*E\p = ip{v),v G A,p = ip{v),v G A, g G A} 

Then Wi C TMi is a constant and therefore integrable distribution whose rank is 
2 dim A. Write for the pullback of the presymplectic form u = dq A dp to 

each integral leaf. Let us compute ker oj^^-^^ . 

Consider (q,v,p,q,v,if{v)),{q,v,p,q',v',ip{v')) G Wi, where q,v,q',v' G A are 
arbitrary. Assume {q,v,p,q,v,ip{v)) G ker a;_^_|_^^ , that is, 

= uj{{q,v,p,q,v,ip{v)), {q,v,p,q ,v' ,ip{v))) = q{ip{v)) - q{^{v)) 

for any q',v' G A. This means that q,v £ {ip{A))°, but since they also belong 
to A then q,v e {<f{A))° D A = {(p{A) + A°)°. That is, ker a;^_^,;^^^ is defined by 
q,v £ {^{A) + A°)°, p = ip{v). The subspace {(p{A) + A°)° is generated by the 
loop currents on the non-inductive loops (Theorem 6.4). Since there are no purely 
capacitive loops, then it is generated by the loop currents on the empty loops. As 
a conclusion, x + Wi is symplectic if and only if there are no empty loops on the 
circuit. 

If there is at least one purely capacitive loop, then the final constraint subman- 
ifold is M3. Then 

W3 = {iq,v,p,q,v,p) £TEeT*E\{q,v,p,q,v,p) £ TM3 , g G A} 

= {{q,v,p,q,v,p) G TE®T*E\p = ip{v),v G A,g G V"^((/?(A) + A°), 

V G V"Hv'(A) + A°),p = ^{v),v G A,4 G i^'H^iA) + A°), 

V £^-\ipiA) + A°),q£ A}. 

Again, W3 C TM3 is a constant and therefore integrable distribution whose rank 
is 2dim(A n ^"^ (93(A) + A°)) = 2dimAi. Write oj^_^_^y.^ for the pullback of the 
presymplectic form uj = dq A dp to each integral leaf, and let us compute ker . 

Consider {q,v,p, q,v, ip{v)), {q,v,p,q' ,v' , Lp{v')) G W3, where q,v,q',v' G Ai are 
arbitrary. Assume {q,v,p,q,v,(p{v)) G kerw^.^^^, that is, 

= uj{{q,v,p,q,v,ip{v)), {q,v,p,q' ,v' ,ip{v'))) = q{(p{v')) - q'{^{v)) 

for any q' ^v' G Ai. This means that q,v £ ((/9(Ai))°. If there is an empty loop, 
represented by 77 7^ 0, then r/ G A, ip{7]) = and r/ G Ai. Also, (p{ri) = 0, so 
{q, v,p, rj, T], 0) is a nonzero element of kera;_^_|_^^(g, v,p) and x+M/3 is not symplectic. 
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If there are no empty loops, we apply the following argument to g, and the 
same conclusion holds for v. Since q E {ip{Ai))° n Ai, then q{(p{q)) = 0, which 
means that the components of q corresponding to the inductive branches are zero, 
as we reasoned in the proof of Theorem 6.4. Also, G Ai C A, so G Ani, 
that is, it can be regarded as a branch current assignment on the non-inductive 
subcircuit, satisfying KCL. In addition, (ji/Ci satisfy KVLNI, which means that 
(ji/Ci satisfy KVL for the non-inductive subcircuit. Now we will apply the fact that 
the dynamics of a circuit with no inductors and no empty loops consists only of 
equilibrium points. Indeed, KVL for such a circuit is a homogeneous linear system 
on the charges qi of the capacitors, so it has a solution q = {qi, . . . ,qn)- Setting 
q{t) = q and v{t) = q(t) = 0, we have that v satisfies KCL trivially. If there was a 
nonzero loop current r], then it would produce a change q = r] in the charges of the 
capacitor on that loop. Then they would no longer satisfy KVL for that loop, since 
d/dt{r]{tp{q))) = i]{ip{r])) ^ 0. Then the currents must therefore be zero, that is, 
q = 0. Of course, this is not true if there are empty loops, which can have arbitrary 
loop currents. By the same reasoning, i; = 0, so kei uj^_^_yy^ = and x + W3 is 
symplectic. ■ 

Let us now write equations of motion in Hamiltonian form, by working on a 
symplectic leaf. This means that the initial conditions must belong to that particular 
leaf. This is related to the comments on equations (5.1) and (5.9) made before. An 
equation of motion in Poisson form that is valid for all initial conditions in Mc will 
be given in Theorem 7.4 by the equations (7.10) and (7.11). 

First, let us represent the KCL equations hy ki, Ka £ E* , in the sense that 
Atj(t') = represents the KCL equation corresponding to node i. The number a 
is the number of independent nodes (usually the number of nodes minus one). As 
before, denote the KVLPC equations by rji, . . . ,r]f, G E, where b is the number of 
independent purely capacitive loops. Seeing each Ki as a row vector and each r]j as 
a column vector, write 



K 



r 



The space TE ®T* E with the presymplectic form a; is embedded in T*TE with 
the canonical symplectic form (see Appendix A), where the variables are {q, v,p, v). 
The integral leaf xq + Wc, where c = 1 or c = 3 as before, and xq = {qo,vo,po), is 
defined as a subspace of T*TE regularly by the equations 



i^i{q) - Ki{qo) 
Pi - LiVi 

Ki{v) 
Vi 



0, 
0, 
0, 
0, 
0, 
0, 



,a 

,n 
,a 

,& 
,n 
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and the matrix of Poisson brackets of these constraints is, in block form, 



S = 









K 









































-if 

















K 
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where (p and ip stand for the diagonal matrices with Li and 1/Cj along their diago- 
nals, respectively. 

A concrete example. We shall illustrate our method with the simple LC circuit 
studied in Yoshimura and Marsden [2006a] which is shown in Figure 6.1. It is a 
4-port LC circuit where the configuration space is E = M.'^. 

f nmp T 1 





1 ^2 1 


1 — ^ — 1 






A 

VC2 '"Ci ^ 











Figure 6.1: The constraint algorithm is apphed to this 4-port LC circuit. 

We shah use the notation q = {qL,qc-L,QC2^<lC3) ^ E, v = {vl,vc^,vc2,vc-s) ^ 
TqE, and p = {pl,PCi,PC2iPC3) S T*E. The Lagrangian of the LC circuit is 
C: TE R, 

.2 Hqc,? Hqc? Hqc^? 

C{q,v) - -L{v,) -2^-2^-2^- 

We will assume the physically meaningful case where < L and < Cj < 00, 
i = 1,2,3. Then, by Theorem 6.6, the CAD algorithm will stop at most at M3. 
Also, Theorem 6.4 gives a direct description of the specific equations that define 
Ml, M2 and M3. This circuit has a purely capacitive loop, so Mi D M2 2 -^3- 
However, the algorithm could also be applied for the case where negative values for 
L or Ci are allowed. For this particular circuit, the algorithm stops at M3 regardless 
of the signs of the inductance and capacitances. 

The KCL constraints A C TE for the current v are 



(6.29) 
(6.30) 
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Therefore the constraint KCL space is defined, for each q G E, hy 
Aiq) = {v£ T,E I (p'^^v) = 0, a = 1, 2}, 

where 

^Pk) - 1^ -11-1 

On the other hand, the annihilator A° of A is the constraint KVL space, defined, 
for each q £ E,hy 

A°{q) = {pG T*E I (p,v) = 0, for all v G A{q)}. 

Note that {p", a = 1,2} is a basis of A°. Taking into account that the Dirac 
structure Da ^ TM®T*M on M = TE®T*E associated to the space A{q) is, for 
each {q,v,p) £ M, given by 

DA{q,v,p) 

= {{q, v,p, q, v,p, a, 7, /3) G TM T*M \ q E A{q), a+pG A°(g), /3 = 7 = 0}, 

and the corresponding energy <S: M — >■ M, is given by £{q,v,p) = pv — C{q,v), we 
can easily verify that the Dirac dynamical system 

{q,v,p,q,v,p) e d£{q,v,p) € DA{q,v,p) 

is equivalent to the IDE system that is given in coordinates by (6.20)-(6.23), which 
gives 

qi = VL, qci = vc-,, qc2 = VC2, qcs = VC3 (6-31) 
PL + ^ +m + 7^ +m = (6.32) 

PL + ^+PC2 + ^+PCs = (6.33) 

PL = LvL, PCi = PC2 = m = (6-34) 
VL = VC2, vci = - VC3- (6.35) 

We now apply the constraint algorithm CAD for Dirac dynamical systems. We 
calculate the expressions of 

V'(^) = ^ and v(g) = 1^ 



to get 



V'{vL,vci,vc2,vcs) = (LuL, 0,0,0) 
i^{qL,qc^,qc2,qc3) = (0,-7^,-7^,-7^, 

Oi O2 O3 
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From now on we will consider the constant distribution A as a subspace of E. As 
we have seen before Mi = {{q,v,p) \ v £ A, p — ip{v) = 0}, so using (6.29)-(6.30) 
and the expression for ip we get 

Ml = {{qL,qci,qc2,QC3,VL,vci,vc2,vc3,PL,PCi,PC2,PC3) I 

PL = LVL, PCi = PC2 = PCs =0,VL = VC2, Vc^ = VC2 - VCs}- 

We could calculate M2 using the expression (6.24). However, according to Theo- 
rem 6.4, M2 adds one more constraint, corresponding to the KVL for the purely 
capacitive loop that involves capacitors Ci and C3. Then 



PL = LvL, 
QC^ _ QC3 



M2 = \^{qL,qci,qc2,QC3,VL,vci,vc2,vc3,PL,PCi,PC2,PC3) 

PCi = PC2 = PCs =0,VL = VC2, Vc^ = VC2 - vcs, - 

As explained right after the proof of Theorem 6.4, the final constraint submanifold 
M3 adds the constraint /Ci = vcs / C3 : 

M3 = \^{qL,qCi,qC2,QCs,VL,vci,vc2,vcs,PL,PCi,PC2,PCs) 

PL = LvL, PCi = PC2 = PCs =0,VL = VC2 , = - vcs , 
qci _ qcs vci _ vcs \ 

To solve the system we can simply parametrize M3, which has dimension 4, for 
instance by taking some appropriate 4 of the 12 variables 

ilL , qci , qC2 , QCs , VL , vci , VC2 , vcs , Pl , PCi , PC2 , PCs ) 

as being independent parameters and then replace in equations (6.31)-(6.33) to 
obtain an ODE equivalent to equations of motion. For instance, if we choose qL, 
QCi, qc2i PL 1 sts independent variables we get the ODE 



QL 



PL 

L 

Ci 



qc2 

PL 



PL 

L 

Qci _ qch 
Ci C2 



We know that solutions should be tangent to VFs = TM^nE^^ where, according 
to (6.10), E^^ is defined by the conditions q £ A, that is, 

<iL = qc2, QCi = qc2 - QCs- (6-37) 

Therefore W3 is defined by (6.37) together with the equations defining TM3 which 
are obtained by differentiating with respect to time the equations (6.36) defining 
M3, that is, 

r . . . . n QCi qCs vc, VC3 

PL = LvL, PCi = PC2 = PCs = 0, VL = VC2, vc^ = VC2-VCS, -PT = -PT^ 7^ = 7^" 
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Then W3 is a constant and therefore integrable distribution of rank 2. Its integral 
leaves are preserved by the flow. One can check that the pullback of the presymplec- 
tic form lo = dq A dp to each leaf is nondegenerate, so they are symplectic manifolds. 
Denote W3 = W^i^O), so these leaves can be written as x + W3. 
For a given xq S M3, say 

2^0 = iQL0,QCi0,QC20,QC30,VL0,Vc\0,VC20,VC30,PL0,PCi0,PC207PC30) 

we have that elements of the symplectic leaf xq + W3 are characterized by the 
conditions defining M3 plus the condition obtained by integrating with respect to 
time the conditions (6.29) and (6.30), applied to {qL—QL0,QCi~QCi0,QC2~'lC20yQC3 — 
qcso), that is 

-{ql - QLo) + {qc2 - QC20) = 

-{qci - qcio) + {qc2 - 1020) - {qc-i - qcso) = 0. 

For simplicity, we will assume from now on that the conditions 

qw - qc20 = 
qcio - qc20 + qcso = 0, 

are satisfied. While this is not the most general case, it is enough to illustrate the 
procedure. 

We can conclude that xq + W3 has dimension 2 and in fact the projection 7f(xo + 
W3) C T*E is simply the subspace of T*E defined by the conditions qci = qc2 = 
qc^ = and Pd = PC2 = PC3 = 0. Therefore one can use the variables {qL,PL) to 
parametrize 7f(xo + W3), and, in fact, we obtain 

'"L = qci = (-j_^^fj^ QL (6.38) 

vci = (j^_^f-,J lPL qc2 = qL (6.39) 



PL C3 
VC2 = Y = (j^^(j^ lL (6.40) 



^C-s = TT-TT^TPi PCi = PC2 =PC3=0 (6.41) 



= z^Ci = t^Ci = T^Cz = 0. (6.42) 
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Let xq + W3 be defined by ej = 0, j = 1, . . . , 14, where we have chosen 





= PL- 


- LVL 
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= Dr" 




£3 


= PC2 






= PC3 




es 


= VL- 


- VC2 


ee 


= vcj 


- VC2 + 




qci 


QC3 




Ci 


C3 



We can calculate the matrix (Sj 





_ J_ 

Ci 


_!_ 
C3 






















which is invertible, with inverse S*^ then S'^^'Sfej = 6). This means that xq + W3 is 
a second class submanifold, or equivalently, that all the constraints are second class. 
In particular, xq + W3 is a symplectic submanifold, and we will call the symplectic 
form Uc- Then we have two ways of writing equations of motion on this submanifold. 

On one hand, we can write the equations of motion on the symplectic manifold 
Xq + W3 in the form (2.3), 
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{ujcf {d{£\xo + W^)) . 



■^S\xo+W-i 

For this we are going to use the previous parametrization of X0+W3 with coordinates 
qL, PL, so = dqi /\dpL- Since the Lagrangian is given by 



C{q,v) 



1 



i(gcj' Hqc2? HqcJ' 



2"'^^' 2 Ci 2 C2 2 C3 ' 
in terms of the parametrization the energy £ = pv — C{q, v) restricted to xq + W3 is 



1 



1 



+ 
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Then the vector field X^^^^_^y^_^ is given in coordinates qi, Pl by 

with H = 8\xq + W^, for short. 

We remark that the previous equations can easily be deduced by elementary 
rules of circuit theory. Namely, one can use the formulas for capacitors in parallel 
and series and replace the three capacitors by a single one. The resulting circuit is 
very easy to solve, and the corresponding system is equivalent to equations (6.43) 
and (6.44). However, we must remark that this simplified system no longer accounts 
for the currents and voltages on the original capacitors. 

On the other hand, we can find the evolution of all the variables {q,v,p) using 
the extended energy (3.26) rather than the total energy, since xq + W3 is symplectic 
and using the last paragraph of remark 3.30. Namely, we will calculate the vector 
field X = d£^ + A-^g^Xf^. associated to the extended energy. 

We have A^^g^ = T,^^{£, ej}, j = 1, . . . , 14, where the column vector [{£, ej}] is 



[{£,ej}] 



0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 



Oi O2 C/3 

For instance, we can calculate pL and vl as 

PL = X{pl) = {d£f{PL) + AJ3)X,^.(PL) = {PL,£] + X\^){PL, e,} = + A93)(-1) 

_ 9Ci qc2 



Ci c, 



2 

VL = X{vl) = {d£f{vL) + \{^)X,^{vl) = {VL,£} + Aj3){^;L,e,} = + 

= 1 _ ^ 

L\ Ci C2 
The complete coordinate expression of X is 

X = {qLACi, qC2 , QC3 , VL , VCi , VC2 , VC'i , PL , PCi , PC2 : PCs , Z>L , i>Ci , %2 . %3 ) 

1 f qci qc2\ Ci ^ ( qci qc2 



1 f qc. qC2 \ C, I ( qc, qc2 \ qc, '^^^ 0, 0, 0, 0, 0, 0, 



L \ Ci C2 J Ci + C3 L \ Ci C2 J Ci 

Note that this is not a Hamiltonian vector field. However, if we use Theorem 3.28, 
we can write the equations of motion in terms of the Dirac bracket and the abridged 
total energy corresponding to the chosen leaf xq + W3. In section 7 we will return 
to LC circuits and show how to write equations of motion on Mc without specifying 
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any leaf a priori, by using an abridged total energy and the Dirac bracket associated 
to adapted constraints. 

We remark that it is not always the case that there is uniqueness of solution. 
In fact, several branches of the circuit may be present with no capacitors or induc- 
tors, as a limit case. In the next section we show how to write a Poisson bracket 
description of the solution. 

7 An extension of the Dirac theory of constraints 

In sections 2 and 3, under precise regularity conditions, we studied several geometric 
and algebraic objects related to the Dirac theory of constraints like primary and final 
constraints, first class and second class constraints, Dirac brackets, first class and 
second class constraint submanifolds. 

The notions of first class and second class depend only on the final constraint 
submanifold Mc and they do not depend on the primary constraint submanifold or 
the Hamiltonian, which are needed only to write equations of motion. Then, as we 
have already said at the beginning of section 3.3, in order to study those notions one 
can start with an abstract situation, given by a symplectic manifold (P, Q) and an 
arbitrary submanifold 5 C P, to be thought of as the final constraint submanifold. 

In this section, motivated by integrable nonholonomic mechanics and LC circuits 
theory, we extend those studies by replacing 5 by a submanifold S C P regularly 
foliated by submanifolds S. This means that S = UsgeS" and the map S — >• © given 
by X I—)- 5 iff x G is a submersion. To study the dynamics, one should consider 
a primary regularly foliated constraint submanifold S' whose leaves S' satisfy the 
condition S = 5" H S, and an energy 8: P — )• M. 

The meaning of S can be interpreted in connection with the constraint algorithm 
CAD for a Dirac dynamical system (1.1) on the manifold M in the case in which the 
Dirac structure D is integrable, which gives a foliation of M. We should compare 
S with the foliated submanifold Mc appearing in remark 5.3, while S' should be 
compared with M submanifold of P. 

The results in section 3.4 can be extended for a certain kind of Dirac dynamical 
systems. Namely, let S C S' C P as before. Consider each leaf S'^ endowed with the 
presymplectic form ujc obtained as the pull-back of Vl. Assume that the distribution 
ker ujc is regular and its dimension does not depend on C. Define the Dirac structure 
D on S' by declaring that its presymplectic leaves are {S'^^uc)- We will consider 
the Dirac dynamical system (1.1) for the case M = S', that is, 



The results in section 3.4 would correspond to the case in which the foliation of S 
consists of a single leaf. 
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7.1 Dirac brackets adapted to foliated constraint submanifolds 

Description of primary and final foliated constraint submanifolds. Let 

{P, be a symplectic manifold. We are going to describe the primary foliated 
constraint submanifold S' and the final foliated constraint suhmanifold S of P. We 
will assume that there is a one to one correspondence between the foliations such 
that for each leaf S' C S' the corresponding leaf of S is 5" = 5" n S. A motivation 
for this assumption comes from mechanics, where an energy on P is given. In 
that context, each leaf 5 C S is the final constraint submanifold corresponding to 
the primary constraint submanifold S' C S' after applying a constraint algorithm. 
These leaves will be parametrized by a certain vector C, and we will denote them 
by S'(j and Sc respectively. We are going to work on a suitable neighborhood U of 
S' in P, where all these submanifolds can be defined regularly by equations, as we 
will describe next. 

For convenience, we will work with a more general situation where S' is a member 
of a family of primary foliated constraint submanifolds parameterized by a certain 
vector C(i,a') ^iid S is a member of a corresponding family of final foliated constraint 
submanifolds parameterized by a certain vector C(i q). For C(i = and C(i = 
we obtain S' and S, respectively. 

Let 

C = (Ci, . . . , Ca'.Ca'+l-, . ■ • , Cq, Ca+1, . . . , Cfc) = {C (i^a')-, C l^a' +l,a)^ Ci^a+l,b)) (7-2) 

be a generic point of M'' = W' x W'"' x M^"", where a' <a<h. Denote C(i^a) = 
(Ci, . . . , Co)- Then we assume the following definitions by equations (regularly on 
U) of the submanifolds S*^(i''»') and S*"^'!-") and their foliations by leaves 5^ and Sc 
respectively, which satisfy Sc = S'^j r\ S'^t^-") , namely, 

= Cj, i = 1, . . . , a', defines S^d-"') (7.3) 

(/,. = Ci, i = 1, . . . , a', a' + 1, . . . , a, defines 8*^(1.-) (7.4) 

4>i = Ci, i = 1, . . . , a', a + 1, . . . , 6, defines S'c (7.5) 

(pi = Ci, i = 1, . . . ,a,a + 1, . . . ,b, defines Sc (7.6) 

Since for each primary constraint submanifold S'q there is only one final constraint 
submanifold Sc, then Ci^a'+i,a) must be a function of C(i and without loss of 
generality we will assume that this function maps to 0. Then for fixed C(i,a')) 
C(a'+i,a) is also fixed and one has, for each choice of C(^a+i,b) = (C'a+i, • • • , Ch), one 
leaf S'^, of S^(i>"') and the corresponding leaf Sc = S'^, n S^d--) of S^d-"). 

Assumption 7.1. From now we assume that our foliated submanifolds S' and S 
are S t^'"'' and S (i-") for the choice C(i = and C(i = 0, respectively, while 
the leaves S' and S are the leaves S'f-, and Sc respectively, for the choice C = 
{0,0,C(^a+i,b)) ! o.'f^d C(q_|_i (,) varying arbitrarily, that is, C G {0} x {0} x M*"" (see 
Figure 7.1). 



7.1 Dirac brackets adapted to foliated constraint submanifolds 64 



S' 




Figure 7.1: S' is foliated by the submanifolds {S'(j}, while S is foliated by {Sc}- Here we 
depict the leaves for two values C and C. 

Remark 7.2. The case considered by Dirac, studied in section 3, corresponds to 
a = b. It is in this sense that our theory extends Dirac's. In this case, the primary 
constraint S' and the final constraint S have only one leaf. More precisely, this 
particular case has been studied in Lemma 3.23, with a different notation. Our 
extension takes care of those Dirac dynamical systems coming for instance from 
integrable nonholonomic systems or LC circuit theory. 

Equations of motion as a collection of Hamilton's equations. Let us adopt 
the point of view that the Dirac dynamical system (7.1) is a collection of Gotay- 
Nester systems, one for each leaf, namely 

n{x){x,6x) = d£{x){6x), (7.7) 

where {x, x) G T^Sc, for all 5x G T^S'^j. 

Then one can apply the theory developed in section 3.4 leaf by leaf. For each 
C = (0, 0, C(a_|_i b)), the total Hamiltonian defined in (3.19) should be replaced by 
the total energy depending on C, by definition, 

£c,T = £ + y[4>i ~ Ci), sum over i = 1, . . . ,a' ,a + 1, . . . ,b 

and, for each such C, the A* G C°°(p(), i = 1, . . . , a', a + 1, . . . , 6 must satisfy the 
conditions {£c,T-, (pj ~ Cj}{x) = 0, j = 1, . . . ,b, for each x £ Sc or, equivalently, 

{£, (Pj}{x) + X'{<Pi, 0,}(x) = 0, j = 1, . . . , 6, 

for each x £ Sc- This shows in particular the interesting fact that the set of 
solutions, which is an affine space, depends on x but not on C. We assume, as 
part of our regularity conditions, that the solutions (A^, . . . , A" , A""*"^, . . . , A*) form 
a nonempty affine bundle A — >• S. Since we are working locally, by shrinking U 
conveniently we can assume that this bundle is trivial. Furthermore, we extend this 
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bundle arbitrarily to a trivial affine bundle over which we continue to denote 
h. ^U. Note that for each C, the total energy can be seen as a function £c,T ■ A — )• 
M.. ^ 

For each section A of A one can define Sq^j^ = X*£c,T : W — 5- M, that is, 

£^jj,{x) = £{x) + X^{x){(pi{x) — Ci). One has an equation of motion on U, 

X£^--^^[x) = [d£^~^,^)\x), (7.8) 

which should be interpreted properly, as follows. In the definition of <?c,T; x and C 
are independent variables. For a given xq G let Co = [4>i{xq), . . . , ipbixo)); then 
^Co A t(^) ^ function of x and 

^^co,a,t(^o) = (^^Co,a,t)"(^o) = (xo) + AHxo)#f (xo) . 

This means that equation (7.8) should not be naively interpreted as Hamilton's 
equation, and therefore the algorithms devised by Dirac and Gotay-Nester cannot 
be generalized in a direct way for the foliated case. In fact, even though the variable 
C is a function of x, it should be considered a constant in order to calculate the vector 
field at the point x. Note that, by construction, the condition Xs - (xq) G Tx^Scn 
is satisfied, which means that the leaves Scq of S are preserved by the motion, as 
expected. 

Remark 7.3. For any given solution x{t) of equations of motion (7.7) there exist 
uniquely determined AJ,z = l,...,a',a + l,...,6 such that 

x{t) = d£^ {x{t)) + Xid4 {x{t)) . 

On the other hand, for any given time-dependent section AJ(x), i = 1, . . . ,a',a + 
1, . . . ,b of A one has a time-dependent vector field 

d£^{x) + Xi{x)d4{x) (7.9) 

whose integral curves are solutions to the equations of motion. One can show that 
under the strong regularity conditions that we assume in this paper, all the solutions 
of the equations of motion can be represented in this way, at least locally. 

Summarizing, the vector fields (7.9) obtained from (7.8) are not in general Hamil- 
tonian vector fields with respect to the canonical bracket, as we have indicated be- 
fore. Now we will see how to recover the Hamiltonian character of the equation 
of motion. We are going to write it as Hamilton's equation in terms of a certain 
Poisson bracket whose symplectic leaves are isomorphic to U. 

Equations of motion in terms of a Poisson bracket. Let us consider an 
extended Poisson manifold, namely the product Poisson manifold (U x {0} x {0} x 
M^~", { , }(^)) with the canonical Poisson bracket (f^)~^ on U and the Poisson 
bracket on {0} x {0} x M''"". Let £:t: A x {0} x {0} x M^"" M be defined by 
£t{X,C) := £c,TiX). Given a section A of A one can define an energy function 
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£j^j,:U X {0} X {0} X R''-" ^ M as £^^{x,C) := £c,t{Hx)) = £^^j,{x). Consider 
its corresponding Hamiltonian vector field 

onUx {0} X {0} X M^~'^ with respect to the Poisson bracket { , }^''\ where Xs^ - ^ (x) 
is defined in (7.8). We remark that even though the first component cannot be 
interpreted in general as a Hamiltonian vector field on U, as we have observed 
before, Xf- is Hamiltonian on the extended Poisson manifold. 
Consider the function 

C(o,o)(^) = (0, 0, </.a+i(x), . . . , Mx)) G K"' X M"-"' X M^-" 

defined on U. Clearly, graph(C(o,o)) is a submanifold oi U x {0} x {0} x R^^"- 
diffeomorphic to lA. 

Since X^^-^ ^ preserves Sc for each C, it is straightforward to prove that X^^^ ^ 
preserves graph(C(o,o)l'S'c)5 which for each C is a copy of the level set C^~^q^(C) 
embedded in ^ x {0} x {0} x M''^'^, and therefore X^^^ preserves graph(C(o,o) 
which is a copy of S. 

Consider the evolution equation for a quantity F onU x {0} x {0} x M^~" given 

by 

Fix,C) = {F,£~^^^}^'\x,C). 

Note that this is an equation in Poisson form on U x {0} x {0} x M^^", not on 
S X {0} X {0} X M}-"'. Note that for a = 6 and P = T*Q, then S = 5 and we recover 
the situation in Dirac's theory, where Hamilton's equations of motion are written 
in a neighborhood of the final constraint submanifold in T*Q rather than on the 
constraint submanifold itself. 

Local equations of motion in Poisson form with respect to the Dirac 
bracket. As we have remarked above, the equations of motion in terms of the 
canonical bracket for the foliated case are not naive extensions of the Dirac or 
Gotay-Nester procedures. A similar situation occurs with the equation g ~ [g, Ht]* 
in Dirac [1964], page 42. As we will show next, we can extend this equation to the 
foliated case by using adapted constraints and the abridged total energy. 

In the definition of S'q and 5c take C = 0. Then we can apply the pro- 
cedure of section 3.4 with S = Sq and S' = S'q, and obtain second class con- 
straints xi, • • ■ ) X2s among (pi, i = 1, . . . , 6 adapted to Sq and also we can choose 
some primary constraints among 0j, i = 1, . . . , a', a -|- 1, . . . , 6, say w.l.o.g. (/>fc, k = 
1 , . . . , a' — , a -f- 1 , . . . , 6 — s'^ . Observe that here s' = s'j^ + s'^ has the same meaning 
as in section 3.4, this time considering the constraints 0i = 0, ...,(/>(, = 0. Then the 
equations of motion on Sq can be written in the form (3.25). 

Now, for any small enough C, which defines Sc, one can readily see that xi ~ 
Bi, . ■ ■ , X2s — are second class constraints adapted to 5^ where i?i , . . . , i?2s 
are the components of C corresponding to the choice xii ■ ■ ■ iX2s- This implies 
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immediately that the Dirac bracket for Sc does not depend on C. On the other 
hand, each (j)/^ — C^, k = 1, . . . ,a' — s'^, a + 1, . . . , 6 — s'^, is a primary constraint for 

S'c- 

Using the previous facts and the fact that the A''^ are arbitrary (even time- 
dependent) parameters, as it happens with the A'* in equation (3.25), one can con- 
clude that the equation of motion on each Sc is given by 

F = {F,£}^^) + X''{F,cl>k}i^) 

= {F,£: + A'Va(x)> (7-10) 

sum over k = 1, . . . ,a' — s'^, a + 1, . . . ,b — s'^. 

This means in particular that each Sc is preserved by the motion, or that the 
Hamiltonian vector fields 

-^(x),£-+A'fe</,fe> (7-11) 
defined on U, are tangent to the Sc- 

Note that £ + X'^cpk is the abridged total energy Sat for Sq and S'^ as defined 
in (3.23). 

In conclusion, using adapted constraints, the abridged total energy and the Dirac 
bracket yields a very simple procedure for writing the equations of motion, because 
one only needs to write the abridged total energy for one leaf, say C = 0, and (7.10) 
gives the equation of motion for all nearby leaves. 

We have proven the following theorem, which extends Theorem 3.28. 

Theorem 7.4. Let (P, O) be a symplectic manifold and let S' 3 S 6e given pri- 
mary and final foliated constraint submanifolds, with leaves S'q and Sc = S'^-, (1 S, 
respectively, as described by (7.2)-(7.6). Assume that all the hypotheses of The- 
orem 3.28 are satisfied for the submanifolds S'q 5 5*0. In addition, suppose that 
the number s' appearing in Assumption 3.26 is the same for all Sc for C close 
enough to 0. From the validity of Assumption 3.27 for S'q, it is immediate to see 
that the number 2s is also the same for all S'q for C close enough to 0. Choose 
second class constraints (x'l, ■ ■ . ,x's'^ x's'+i^ • ■ • > X2s) (adapted to S'q and also (pk, k = 
1, . . . ,a' — s'^, a + 1, . . . , 6 — s'^, as in Theorem 3.28. 

Let iS: P — ?• M be an energy function, and consider the Dirac dynamical sys- 
tem (7.1). Define the abridged total energy Sat = S + X'^cj)^ for Sq and S'q as 
in Theorem 3.28. Then, each € S has an open neighborhood U such that the 
vector field (equation (7.11)) defined on U, when restricted to the final 

foliated constraint submanifold S CiU, represents equations of motion of the Dirac 
dynamical system on SOU. Moreover, the evolution of the system preserves the 
leaves Sc O U. In addition, equation (7.10) gives equations of motion on each 
Sc n U . Since X(^) is tangent to Sc n U , then the evolution of a function 
f on Sc U is given by (7.10) for any F such that f = F\Sc n U . Also, each 
Sc U is contained in a unique symplectic leaf of the Dirac bracket, defined by 

x'l = Bl,-- - , X's' = Bs',Xs'+i = Bs>^i, x'2s = B2s- 

Remark 7.5. Dirac [1964, page 42], wrote the equation 

g^{g,HTr 
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which represents the dynamics on the final constraint submanifold, in terms of the 
total Hamiltonian and the Dirac bracket. This equation does not gives the correct 
dynamics for the case of foliated constraint submanifolds. As we have seen, for the 
foliated case, one has, instead a similar equations of motion (7.10) in terms of the 
Dirac bracket and the abridged total energy Sat- 



The concrete example of an LC circuit of section 6 revisited. Consider 
the symplectic leaf S of Mj, defined parametrically by equations (6.38)-(6.42). The 
parameter space carries the canonical symplectic form dqi A dpL, which coincides 
with the pullback of the canonical symplectic form on T*TQ to the parameter space. 
Since the matrix S is invertible, all constraints are second class constraints. The 
number b — s'^ in equation (7.10) is 0, then the abridged total energy is simply the 
energy 8. 

Then it is easy to calculate the Dirac bracket as the canonical bracket on the 
parameter space of the functions (6.38)-(6.42). Let y = (q,v,p,i'), then the matrix 
of Dirac brackets of these variables can be written in block form as 



04x4 A 04x7 
—A'^ OsxS 05x7 
07x4 07x5 07x7 



where 



1 


Ci 


1 


C3 


1 


L 


(Cl+C3)i 


L 


{Ci+C3)L 


Ci 


q2 


Ci 


C1C3 




{Ci+C3)L 




{C\+C3)L 


{C1+C3YL 


(C1+C3) 


1 


Ci 


1 


C3 


1 


L 


(Ci+C3)L 


L 


(Ci+C3)L 


C'i 


C1C3 


C3 


cl 


Cs 


(Ci+C3)L 




{Ci+C3)L 


(Ci+C3)2L 


(C1+C3) 



The evolution equations for the vector 
(2/\ • • • ,y^^) = {qL,qci,qc2,QC3,VL,vcj,vc2,vc3,PL,PCi,PC2,PC3,^L,^Ci,J^C2,i'C3) 



with initial condition belonging to M3, that is, satisfying (6.36): pL = Lvl, pci = 
PC2 = PC3 = 0,VL = VC2, vci = VC2 - VC3, "cf = cf' cf = cf ' IS gi^^^ t>y the 



product of the matrix {{y\ y-'lx) by the vector 
dS 



■ ■ ■ ' 
(o, 



' Ci 



QC2 
'C2' 



qc3 
C3 



PL - VL,PCi,PC2,PC3,VL,vci,vc2,vc3,0,0,0,0 



8 Conclusions and future work 

We have shown that Dirac's work on constrained systems can be extended for cases 
where the primary constraint has a given foliation. This extends the applicability of 
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the theory for Dirac dynamical systems, hke LC circuits where the constraints may 
come from singularities of the Lagrangian and, besides, from Kirchhoff 's Current 
Law. Throughout the paper we combine ideas of an algebraic character from Dirac's 
theory with some more geometrically inspired ideas from Gotay-Nester work. In 
particular, we use Dirac brackets adapted to the foliation as well as a Constraint 
Algorithm for Dirac dynamical systems (CAD), which is an extension of the Gotay- 
Nester algorithm. The results were proven locally and under regularity conditions. 
It is our purpose to study in the future the globalization of the results of the paper 
as well as the singular cases, using IDE techniques, with applications. 
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A Appendix 

Lemma A.l. There is a canonical inclusion ip: TQ © T*Q — t- T*TQ. In addition, 
consider the canonical two-forms ujt*q and ujt*tq on T*Q and T*TQ respectively, 
the canonical projection pi rptQ-. TQ®T*Q — )• T*Q, and define the presymplectic two- 
form Lo = pr^*Q ujt*q on TQ®T*Q. Then the inclusion preserves the corresponding 
two-forms, that is, oj = ip*coT*TQ- 

Proof. If TQ : TQ — > Q and ttq : TTQ — > TQ are the tangent projections, we can 
consider the dual tangent rhombic 

TTQ 




TQ TQ 




Q 

Define ip:TQ® T*Q T*TQ by (p{vq © a^) G T*^TQ, 

ip{Vq © ag) ■ Wy^ = aq ■ TTQ{Wy^), 

for G Ty^TQ. Here Vq © aq denotes an element in the Pontryagin bundle over 
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the point q E Q. Note that the fohowing diagram commutes. 

TQ T*Q T*TQ 

TQ 



TQ 



Let us see that ip is an injcctive vector bundle map from the bundle pr^^g : TQ © 
T*Q TQ to the cotangent bundle tttq : T*TQ TQ, over the identity of TQ. 
The last part of this assertion follows from the commutative diagram above. 

First, if (p{vq © aq) = ^{v'^, © a'^,) then both sides are in the same fiber T*JTQ = 



T*, TQ, so Vq 



Also, for all G T*TQ we have 

ip{Vq © aq) ■ Wy^ = ip{Vq © a' ) ■ Wy^ 



or 



an 



Ttq{wv^) = a' -TTQiwug 



Since tq : TQ — > Q is a submersion, wc have aq = a'^ and (p is injective. 
Second, cp is linear on each fiber, since 

ip{Vq®{aq + Xaq))-W^^ = {aq + Xa'g) -TTQ^WyJ = ip{Vq®aq)-Wv^ + \ip{Vq®a'q)-Wv^. 

For the second part of the lemma, let us recall the definition of the canonical 
one-form on eT*Q G 0^(r*Q). For a, G T*Q, OT-Qiaq) is an element of T*T*Q 
such that for G TaqT*Q, 

OT*Q{aq) ■ Wag = aqiTTTqiWa,)), 

where ttq : T* Q Q is the cotangent bundle projection. 

With a similar notation, the canonical one-form 6t*tq € ^^{T*TQ) is given by 

OT*TQ{av,) ■ Wa,^ = ay^iTTTTQiWa^J). 

Pulling back these forms to the Pontryagin bundle by cp and the projection pirptg : TQC 
T*Q T*Q, we obtain the same one-form. Indeed, for w G T^^^agiTQ ®T*Q), we 
get on one hand 

(prT*Q &T*Q){Vq ®aq)-w = OT^Qioiq) ■ Tptt*q{w) = aq ■ T{ttq o prj,*Q)(u;), 
and on the other hand 

{(p*OT*TQ){Vq ®aq)-W = 9T*TQ{^{Vq © a,)) • Tip{w) = 

^ivq © aq) ■ Ttttq{T(p{w)) = aq ■ T{tq o tttq o (p){w). 
However, the following diagram commutes 

TQ®T*Q 

P^T*Q 

T*Q 
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so ttq o pr rp*Q = tqo tttq o ^ and therefore pr^-.Q 6t*q = ^*Gt*tq- Taking (minus) 
the differential of this identity, we obtain uj = ip*ujT*TQ- ^ 

In local coordinates, if we denote the elements of TQ®T*Q and T*TQ by {q, v, p) 
and {q^v^p^u) respectively, then ip{q,v,p) = {q,v,p,0). 
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